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ABSTRACT 

Four  pinned-base  steel  arches  with  a  96-inch  radius,  143. 8- inch  span,  and  uniform  cross 
section  were  cold-rolled  from  4M13  sections  and  tested  under  various  static  and  dynamic 
loads  uniformly  distributed  over  one-half  the  arc  length.  A  maximum  static  load  of 
72, 000  pounds  was  applied  by  the  NCEL  blast  simulator  using  compressed  air.  A  dynamic 
peak  load  of  64,  000  pounds  was  attained  by  detonating  Primacord  in  the  blast  simulator.  The 
blast  loading  had  a  rise  time  of  about  3  milliseconds  and  a  decay  time  of  about  1.6  seconds. 

An  equivalent  triangular  load-time  function  was  used  for  the  dynamic  analysis.  The  applied 
loads  and  the  resulting  deflections,  strains,  and  reactions  were  measured.  The  reduced  data 
are  presented  in  graphical  and  tabular  form. 

The  theoretical  analyses  for  statically  and  dynamically  loaded  arches  were  based  on  the 
discrete  framework  which  represented  the  continuous  arches  tested.  A  16 -bar  system  was 
used  for  both  static  and  dynamic  response  calculations,  and  a  40-bar  system  was  used  for 
natural  mode  and  frequency  calculations. 

In  the  static  analysis,  the  effects  of  stress  amplification,  misalignment,  and  elastic 
supports  on  the  response  of  the  arch  were  considered.  Due  to  the  strain-hardening 
characteristics  of  the  arch  material,  the  idealized  stress -strain  curve  was  represented  by  a 
tril inear  curve  rather  than  by  the  usual  bilinear  stress-strain  curve. 

A  simplified  dynamic  analysis  gave  results  reasonably  close  to  those  from  more 
rigorous  methods;  the  values  were  on  the  conservative  side.  Nearly  complete  correlation 
between  the  theoretical  and  experimental  results  was  obtained. 

This  work  was  sponsored  by  the  Defense  Atomic  Support  Agency. 
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INTRODUCTION 


The  purpose  of  this  report  is  to  provide  information  on  the  theoretical  and  experimental 
behavior  of  two-hinged  steel  circular  arches  subjected  to  uniformly  distributed  static  and 
dynamic  loads  on  one-  half  the  arch.  Many  theories  and  procedures  have  been  advanced  in 
recent  years  for  determining  the  elastic  and  elastoplastic  response  of  structures  to  uniformly 
distributed  dynamic  loads.  However,  due  to  the  lack  of  large  dynamic  testing  facilities,  there 
have  been  few  direct  comparisons  of  the  theoretical  and  experimental  behavior  of  large-scale 
models. 

The  blast  simulator  at  the  U.  S.  Naval  Civil  Engineering  Laboratory  can  test  nearly 
full-sized  structural  members  under  transient  loads  up  to  185  psi. 

This  report  covers  the  third  and  last  phase  of  tests  on  structural  steel  members;  the 
first  dealt  with  the  static  and  dynamic  behavior  of  portal-frame  knee  connections1  and  the 
second  with  static  and  dynamic  behavior  of  pinned-  and  fixed-base  portal  frames.  2,3  For  the 
third  phase,  four  steel  arches  of  96-inch  radius  and  143. 8-inch  span  were  designed,  fabri¬ 
cated,  ami  tested.  The  primary  objectives  ui  ihe  tesi  series  were  to  obtain  information  on  the 
behavior  of  members  subjected  to  static  and  dynamic  loads  and  to  examine  the  behavior  and 
failure  mode  of  the  members  under  static  and  dynamic  loads.  The  information  obtained  will 
be  used  to  recommend  changes  in  the  present  design  criteria  of  structures  subjected  to  blast 
loads. 

The  work  was  accomplished  under  Work  Unit  Y-P008-08-02-102,  Blast  Resistance  of 
Structural  Shapes,  and  was  sponsored  by  the  Defense  Atomic  Support  Agency  through  the 
Naval  Facilities  Engineering  Command. 

Wherever  possible,  symbols  used  in  this  report  are  those  recommended  by  the  Defense 
Atomic  Support  Agency.  A  complete  list  of  symbols  appearing  in  this  report  precedes  the 
appendixes. 

EXPERIMENTAL  TECHNIQUES  AND  PROCEDURES 
Design  Procedure  and  Description  of  Arch 

Two  criteria  in  determining  the  size  and  strength  of  the  antisymmetrically  loaded  arches 
to  be  tested  were  that  they  should  fit  in  the  blast  simulator  and  that  the  maximum  load  should 
be  limited  to  approximately  50  percent  of  the  capacity  of  the  simulator.  It  was  determined 
that  the  central  angle  of  the  arch  would  be  97  degrees  and  the  radius  to  the  neutral  axis  would 
be  96  inches.  The  width  between  the  skirts  of  the  simulator  is  8  inches.  In  previous  tests,  it 
was  determined  that  the  most  efficient  width  of  the  test  member  should  be  between  7-3/ 4  and 
7-7/8  inches. 

With  these  dimensions  and  loading  conditions,  the  mechanism  method  of  plastic  design 
was  used.  It  should  be  noted  that  the  load  produces  sides  way  and  that  plastic  hinges 
develop  in  the  vicinity  of  the  quarter  points,  b  and  d,  as  shown  in  Figure  1.  The  instantaneous 
center  of  rotation  then  lies  at  the  intersection,  c,  of  straight  lines  drawn  through  the  hinge 
points,  a  and  b,  and  d  and  e;  the  central  angle  Is  designated  ^Q.  For  simplicity,  the  load,  w, 
has  been  replaced  by  two  resultant  forces,  F}  and  F2,  acting  normal  to  chords  ed  and  dk. 

From  Figure  1,  the  following  geometrical  relations  can  be  oWained: 
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External  work  for  very  small  deflections  is  determined  as  follows: 


WE  =  F1  f  +  F2  * 
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or  by  using  Equation  1,  and  employing  identities 


WE  =  2FR  $  sin  -~ 


Internal  work  is  determined  as  follows: 


From  Figure  1,  it  can  be  shown  that 


Wj  =  2Mp  (i//  +  If/'}  a  4Mp^cos  ~ 


Equating  (3)  and  (4)  and  solving 
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Since  is  equal  to  97  degrees 


M  =  0. 0484  R  w 
P 


Using  a  radius  of  06  inches  and  a  design  yield  load  of  75  x  8  or  600  lb/in. ,  the  plastic  moment 
is  calculated  as 


Mp  -  267, 384  in,  -lb 
The  plastic  modulus,  Z,  is  then 


267. 384 
40,000 


6.68  in.3 


By  referring  to  a  plastic  modulus  table,  a  reasonable  section  was  chosen  as 

4M13  lb 
Z  *=  6.1  in.3 


A  section  with  a  smaller  plastic  modulus  was  selected  for  testing  in  the  blast  simulator 
because  the  width  of  the  top  flange  would  have  to  be  increased  to  about  7-3/4  inches.  After 
the  arch  was  rolled  into  shape,  two  steel  strips  1-7/8  inches  by  5/8  inch  were  welded  to  the 
top  flange  of  the  arches  as  shown  in  section  B-B  of  Figure  2.  The  total  area  of  the  cross 
section  was  6. 16  in.2.  For  this  cross  section,  the  elastic  neutral  axis  was  computed  to  be 
2.64  inches  above  the  base  of  the  lower  flange  and  the  moment  of  inertia  was  computed  to  be 
14. 84  In. 4.  The  plastic  neutral  axis  was  found  to  be  3. 58  inches  above  the  base  of  the  lower 
flange  and  the  plastic  modulus  was  7.15  in. 3.  Needle  roller  bearings  were  fastened  to  the 
transverse  web  stiffener  (as  shown  in  Figure  F-l)  to  reduce  friction  as  the  arch  deflected. 

Loading  System 


The  blast  simulator  (Figure  3)  was  used  to  t^st  the  antisymmetrically  loaded  arches. 

The  load  applied  simulated  a  blast  wave  from  a<  .»r  burst  initially  striking  the  surface  of  the 
arch  at  one  of  the  quarter  points  and  also  perpendicular  to  the  longitudinal  axis.  The  response 
of  the  arch,  as  the  blast  wave  engulfed  the  entire  arch,  was  not  taken  into  account  in  this  test 
program.  For  testing  these  arch*1?,  special  extension  plates,  6  feet  in  height  and  8  feet  in 
length,  were  bolted  to  the  bottc;:i  of  the  blast  simulator  and  tied  together  by  three  sets  of 
massive  collar  beams  to  prevent  separation.  Figure  4  shows  one  side  of  the  extension  plates 


removed,  exposing  the  arch  to  be  tested.  Static  loads  were  applied  by  pumping  air  into  the 
simulator  from  a  compressor.  Dynamic  loads  were  produced  by  detonating  Primacord  within 
the  firing  tube  in  the  simulator.  ®  A  blast  loading  was  produced  which  had  a  rise  time  of 


approximately  2  to  3  milliseconds  followed  by  an  experimental  decay  of  controlled  duration 
(Figure  5).  Leakage  of  the  static  or  dynamic  pressure  was  minimized  by  a  neoprene  seal 
placed  on  top  of  the  flange  (Figure  4). 
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igure  2.  Dimensions  of  the  antisymmetrically  loaded  arch. 
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Figure  3.  View  of  the  blast  simulator. 
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Figure  4.  Setting  up  the  arch  in  the  blast  simulator. 
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Figure  5.  Typical  load-time  curve. 


Instrumentation 

Instrumentation  of  the  test  specimens  consisted  of  strain  gages,  pressure  cells, 
deflection  gages,  and  an  accelerometer  located  as  shown  in  Figure  6.  Information  on  the 
measuring  instruments  is  presented  in  Table  1.  The  recording  equipment  is  shown  in  Figure  7. 
In  addition  to  the  transducers  listed  in  Table  1,  deflection  measurements  were  also  obtained 
by  the  use  of  a  rotating  drum  gage,  and  the  reactions  were  measured  with  a  pair  of 
dynamometers  (Figure  8)  at  each  end  of  the  arch.  The  rotating  drum  deflection  gage  revolved 
at  a  constant  speed,  recording  deflection  and  as  the  arch  deflected,  a  pencil  attached  to  the 
end  of  a  taut  piano  wire  inscribed  a  deflection-time  trace  on  the  graph  paper. 

Test  Procedure 

The  general  test  procedure  for  the  arches  was  as  follows:  the  arch  was  vibrated,  a 
static  load  was  applied  (to  failure  for  one  arch  and  to  approximately  60  percent  of  yield  for 
the  other  arches),  and  a  series  of  dynamic  loads  was  applied  to  the  arches  that  were  not 
statically  loaded  to  failure. 

Two  types  of  vibrational  mode  tests  were  conducted.  The  first  was  to  excite  the  lowest 
natural  frequency  of  the  two-hinged  arch  in  the  deflection  mode  and  the  second  was  to  excite 
the  arch  in  the  compression  mode.  For  both  types  of  tests,  the  arch  was  wedged  upward  by 
using  a  piece  of  timber  at  the  desired  location  (at  the  quarter  point  for  the  deflection  mode  and 
at  the  crown  for  the  compression  mode).  The  timber  was  suddenly  knocked  out  with  a  sledge 
hammer.  The  subsequent  vibrations  were  recorded  on  the  oscillograph;  a  typical  oscillogram 
is  shown  in  Figure  9.  A  smooth,  continuous  trace  indicated  that  the  specimen  was  vibrating 
freely  in  the  blast  simulator;  the  natural  period  of  vibration  was  obtained  from  these  records. 
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arches  AH-3  and  AH-4 


Figure  7.  Electronic  amplifier  and  recording  oscillograph. 


vertical  reactions  at  the  loaded  end  of  the  arch. 
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Figure  9.  Free-vibration  test  of  arch  AH-2. 
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Table  1.  Instrumentation  Used  in  Tests 


Measurement 

Description 

Type  or  Model 

Manufacturer 

Load 

200-psi  pressure  cells 

PA208TCA-200-350 

Statham  Instruments, 
Incorporated 

Acceleration 

200-g  accelerometer 

A5A-200-350 

Statham  Instruments, 
Incorporated 

Deflection 

6-,  10-,  and  15-inch  linear 
potentiometer 

108 

Bourns,  Incorporated 

Strains 

SR-4  electrical  resistance 
strain  gages 

A-5-1,  A-R-l 

Baldwin-Lima-Hamilton 

Corporation 

Whittemore  strain  gages 

2 -inch  type 

Baldwin-Lima-Hamilton 

Corporation 

Other 

Recorder  with  36-channel 
oscillograph  amplifiers 

5-119P4-36 

1-113B 

Consolidated 

Electrodynamics 

Corporation 

The  first  arch  was  tested  statically  to  failure;  the  remaining  arches  were  loaded 
statically  to  approximately  60  percent  of  the  yield  load  and  then  tested  dynamically  in  the 
elastic  range.  This  was  done  so  that  the  elastic  static  response  of  the  arch  could  be  compared 
with  the  elastic  dynamic  response.  Typical  dynamic  test  data  recorded  on  the  oscillograph 
are  shown  in  Figures  10  and  11.  The  final  dynamic  load  on  each  arch  was  in  the  inelastic 
range  so  that  information  on  the  dynamic  response  over  a  wider  range  of  pressure  could  be 
obtained. 

After  a  specimen  had  been  tested  completely,  photographs  of  the  overall  arch,  including 
closeups  of  the  entire  specimen,  were  taken  (Appendix  F). 

COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  RESULTS 

Tensile  Coupon  Tests 

Static  Tests.  Standard  tensile  test  coupons,  2-inch  gage  length,  were  fabricated  from 
the  flange  strip,  the  web,  and  the  upper  and  lower  flanges.  Figure  12  identifies  the  locations 
and  symbols  of  static  test  specimens.  It  should  be  noted  that  the  specimens  were  cut  from  the 
mid-spans  of  the  arches  where  very  little  bending  moment  took  place  during  the  tests.  It  was 
assumed  that  the  specimens  were  representative  of  the  arch  material  before  tests.  Typical 
stress-strain  curves  for  each  location  are  shown  in  Figures  13  through  18. 

Table  2  summarizes  the  results  of  static  coupon  tests.  The  average  values  of  yield 
stress  at  0. 2  percent  offset  were  56, 200,  41,700,  55,100  and  60,500  psi  for  the  web,  flange 
strips,  and  upper  and  lower  flanges,  respectively;  the  average  values  of  ultimate  strength 
were  70,000,  84,400,  70,700  and  72,800  psi,  respectively;  the  average  values  of  the  elastic 
modulus  were  30.0  x  10®,  30,9  x  10®,  27.9  x  10®,  and  26.7  x  10®  psi,  respectively.  It  should 
be  noted  that  the  4M13  section  was  less  ductile  than  the  flange  strips.  This  fact  was  reflected 
in  the  2-inch  elongation  values  of  34.3,  22.0,  and  25.8  percent  for  the  flange  strips,  the  web, 
and  flanges,  respectively;  this  was  due  to  the  fact  that  the  4M13  section  had  been  strained 
beyond  the  strain-hardening  point  during  the  cold-roll  process  (Figure  17).  Further  evidence 
is  shown  in  the  stress-strain  curves  of  various  specimens  of  the  arch.  No  characteristics  of 
a  definite  yield  plateau  are  evident  in  the  stress-strain  curves  of  Figures  13,  15,  and  16  for 
the  web  and  upper  and  lower  flanges  of  the  4M13  section. 
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Measurement 
Strain  ( SGI 


Reaction  (R1 ) 


Calibration  for 
1-inch  Trace  Deflection 

3,090  //in./in. 


59.70  kips 


Strain  (SG2) 


Reaction  (R2 


Strain  (SG6 


Strain  (SGI  1  ) 


Deflection  (D) 


Strain  (SGI 2 


970  U  in./in. 


70.40  kips 


1  t460^in./in. 


0.93  in. 


2,940 /lin./in. 


Acceleration  (A 


Strain  (SG4 


Strain  (SGI 3 


Strain  (SGI 4 


Strain  (SG15) 


Load  (Li 


Load  (L2) 


1 ,970  /tin./in. 


1,640 /lin./in. 


jpjHM 

m 

maas 

M 


79.5  Ib/in. 


location  of  specimens 


e 


AH-1 

AH -3 

AF1S3 

AF3S3 

AF1S4 

AF3S4 

AF1S5 

AF3S5 

AW1S3 

AW3S3 

AH-2 

AH-4 

AF2S3 

AF4S3 

AF2S4 

AF4S4 

AF2S5 

AF4S5 

AW2S3 

AW4S3 

identifying  symbols 


orch 


A 

A 


F  2 
W  2 


-flange  material 

narch  AH-2 

—static  specimen 
^  ^ - specimen  ^1 


L 


web  material 


Figure  12.  Locations  of  tensile  test  specimens  for  the  antisymmetrically  loaded  arch. 


Strain,  e  ( in./in.) 


IMil 


iui 


♦The  material  was  deformed  beyond  the  strain-hardening  point  by  the  cold-roll 
process  during  fabrication  of  the  arches. 
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Dynamic  Tests.  To  determine  the  Increase  in  the  yield  point  of  the  steel  under  dynamic 
load,  tensile  specimens  were  fabricated  and  then  tested  with  the  Laboratory's  50, 000-pound 
dynamic  materials  testing  machine.7  The  machine  has  a  maximum  static  capacity  of 
50, 000  pounds  and  head  velocities  up  to  15  in.  /sec.  With  a  booster  facility,  the  static  load 
capacity  can  be  increased  to  80,000  pounds  and  the  head  velocity  can  be  increased  to  30  in,/sec. 
Statically,  the  piston  stroke  is  4  inches;  at  higher  velocities  it  is  0. 15  inch.  For  standard 
tensile  specimens  (0. 2  in.  2  area),  strain  rates  of  1.5  in./in. /sec  car.  be  obtained. 

As  shown  In  Figure  18,  the  specimens  were  taken  from  both  ends  of  the  arch,  where 
bending  moments  were  small  due  to  hinged  supports.  The  specimens  were  assumed  to  be 
representative  of  the  arch  material  before  tests. 

A  total  of  nine  tensile  specimens,  six  'rom  the  web  and  three  from  the  flange,  were 
fabricated  and  tested.  Dimensions  of  the  test  specimens  a  shewn  in  Figure  10.  To  determine 
the  value  of  the  strain  rate  to  be  used  for  the  dynamic  4r*nv,ue  tests,  the  strain  rates  from 
SG-12  (Figure  6)  were  obtained  from  several  dynamic  test.-;  the  average  strain  rate  was 
0. 9  in./ln./sec.  The  strain  rate  was  maximum  at  this  point  because  the  flew.ro).  and 
compressive  strains  were  of  the  same  sign.  Typical  stress-strain  curves,  are  given  in 
Figures  20  and  21  for  the  web  and  lower  flange,  respectively.  It  should  be  noted  that  the 
stress-strain  curve  of  the  lower  flange  does  not  show  a  definite  yield  plateau. 

Table  3  summarizes  the  results  of  the  dynamic  teste.  The  average  values  of  upper  yield 
stress  were  70, 500  psi  and  89, 700  psi  for  the  web  and  lower  flange,  respectively;  the  average 
values  of  lower  yield  stress  were  69, 100  psi  and  66,  500  psi,  respectively;  the  average  ultimate 
strengths  were  80,600  psi  and  82, 500  psi„  respectively.  The  average  increase  in  the  lower 
yield  point  was  23. 0  percent  for  the  web  and  9. 8  percent  for  the  lower  flange.  The  marked 
difference  between  the  two  sections  was  due  to  the  fact  that  tine  flanges  had  been  stressed 
closer  to  ultimate  stress  than  the  web. 

When  structural  steels  are  stressed  above  60  percent  of  the  ultimate  strength,  structural 
properties  of  the  material  are  changed.  To  illustrate  the  effect  of  strain  hardening,  consider 
a  typical  stress-strain  curve  as  shown  in  Figure  22,  K  the  material  were  stressed  to  A  and 
unloaded,  the  stress-strain  diagram  would  return  to  the  strain  axis  in  practically  a  straight 
line,  AB,  parallel  to  the  initial  straight-line  portion  of  the  diagram.  If  the  stress  were  again 
applied,  the  yield  point  would  be  raised  to  A,  the  diagram  following  the  dashed  line  BAC.  The 
ultimate  strength  would  be  practically  unchanged  but  the  ultimate  elongation  would  be  reduced. 
The  effect  of  the  single  loading  to  A  is  to  increase  the  elastic  strength  and  to  decrease  the 
ratio  of  dynamic  to  static  yield  points.  The  higher  the  point  A,  the  smaller  the  percent 
increase  in  dynamic  yiek.  point. 8 

Static  Tests 

Static  tests  rf  four  circular  arches  were  performed  In  the  blast  simulator  as  described 
earlier;  the  archc  ere  loaded  along  one-half  their  arc  length.  For  the  majority  of  the  teats, 
radial  deflections  and  strains  were  measured  at  the  quarter  points  and  the  crown;  horizontal 
and  vertical  reactions  were  measured  as  shown  in  Figure  8.  Only  arch  AH-1  was  tested  to 
failure;  the  other  three  arches  were  loaded  to  a  maximum  load  of  243  lb/in.  During  the  tests 
it  was  observed  that  the  support  cart  (Figure  8)  moved  a  fraction  of  an  inch  to  the  right;  the 
exact  distance  was  undetermined.  This  was  partially  due  to  the  elastic  strain  of  the  reaction 
dynamometers  and  partially  due  to  the  movement  of  the  tiedown  bolts  for  the  support  cart.  By 
proportioning  the  assumed  total  horizontal  movement  at  the  maximum  load,  the  cart  movement 
was  introduced  in  the  theoretical  analysis  as  a  boundary  condition  by  assuming  that  the 
deformation  was  proportional  to  the  horizontal  reaction  load  at  Rj ;  this  factor  was  taken  to  be 
6.25  x  10-8  in.  /lb. 


Strati,  f  (ki! ) 


Strain,  •  (in./in.) 


Figure  21.  Dynamic  stress-strain  curve  of  the  lower  flange  at 
a  strain  rate  of  0. 517  in./in./sec. 


i'l 


Stroin  (in./in.) 

Figure  iJ2.  Strain-hardening  effect  in  structural  steel. 
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Table  3.  Dynamic  Properties  of  Steel  for  Antisymmetrically  Loaded  Arches 


Location 

of 

Specimen 

Arch  No. 

Upper  Yield 
Stress 

(psi) 

Lower  Yield 
Stress 

(psi) 

Ultimate 

Strength 

(psi) 

Elongation 
in  2  Inches 

<%) 

Strain 

Rate 

(in./in./see) 

Increase  in 
Lower  Yield 
Point  (**.) 

Web* 

AH-2 

74, 200 

72, 100 

83,100 

BETH 

■m 

28.3 

AH-2 

69, 700 

79,000 

17.2 

AH -3 

69, 000 

68,250 

80, 500 

KSSfl 

20.8 

AH-3 

EiKJjKOpS 

83, 500 

BIAS 

29.0 

AH-4 

66,750 

78,100 

25.0 

0.746 

21.4 

AH -4 

69, 100 

67, 800 

79,600 

27.0 

0.652 

23.3 

Average 

69, 100 

24.7 

— 

23.0 

Lower* 

mm 

64,200** 

81,000 

22.5 

Bn 

Flange 

69, 100 

84,300 

— 

I 

AH-4 

mm 

66,200 

82,100 

— 

0.  700 

9.4 

Average 

69, 700 

66, 500 

82,500 

22.5 

— 

9.8 

♦  The  material  was  deformed  beyond  the  strain-hardening  point  by  the  cold-roll  process. 
♦♦Yield  point  stress  at  0. 2%  offset. 


Arch  AH’l  was  first  preloaded  to  487  lb/in.  in  increments  of  162  lb/in.  There  were  no 
appreciable  permanent  strains,  deformation,  or  residual  reactions  observed  after  this  test. 
Data  from  this  test  are  not  given  in  this  report.  Arch  AH-1  was  then  reloaded  to  568  lb/ln. 

At  this  point,  it  was  expected  that  certain  sections  of  the  arch  would  yield.  The  pressure  was 
then  increased  in  small  increments  to  ensure  that  yield  conditions  would  be  recorded.  At  a 
load  of  791  lb/in.,  air  leakage  at  the  neoprene  seal  exceeded  the  compressor  capacity  and  a 
further  increase  in  pressure  was  not  possible.  The  reduced  test  data  are  presented  in 
Table  A-l.  The  maximum  deflection  obtained  was  1. 39  inches;  the  permanent  deflection  was 
0.40  inch  after  unloading.  The  stiffness  of  the  arch  in  the  elastic  range  was  902  lb/in.  /in. 
and  was  close  to  the  theoretical  value  of  937  lb/ln./ln.  The  strain  data  of  Table  A-l  were 
transformed  into  moments,  direct  forces,  principal  stresses,  and  shears;  they  are  given  in 
Table  A-2. 

After  repair  of  the  neoprene  seal  and  retesting,  the  load -deflection  curve  rose  slightly 
higher  than  the  point  of  unloading  before  flattening  out  (Figure  23).  This  was  the  result  of 
strain  aging  of  the  steel.9  When  strain  aging  is  permitted,  the  sharp  yield  and  yield 
elongation  properties  of  mild  steel  return  and  increase  the  entire  stress-strain  curve.  The 
arch  collapsed  at  a  load  of  876  lb/in.  with  a  corresponding  deflection  of  2.  57  inches  at  the 
3/4-point.  No  deflection  was  recorded  at  the  1/4-point.  Figure  F-l  shows  the  collapsed  arch. 
There  was  evidence  of  buckling  at  the  lower  flange  of  the  1/4 -point  as  a  result  of  compressive 
strain.  At  the  3/4 -point  the  lower  flange  and  the  web  were  ruptured. 

The  theoretical  load-deflection  values  (calculated  in  Appendix  B)  are  plotted  in  Figure  23 
for  comparison.  Data  for  this  test  are  given  in  Table  A-3.  Comparisons  of  experimental  and 
theoretical  load-strain  curves  at  the  quarter  points  of  the  arch  are  presented  in  Figures  24 
and  25.  Excellent  correlation  with  strain  gage  4  was  obtained.  Other  theoretical  strain 
curves  were  very  close  to  their  respective  experimental  curves  in  the  elastic  range  but 
deviated  somewhat  in  the  inelastic  range.  A  comparison  of  theoretical  and  experimental 
direct  forces  in  the  arch  is  shown  in  Figure  26.  The  experimental  values  showed  that  the 
direct  forces  at  the  l/4-point  and  at  the  crown  were  very  close  to  each  other  but  were  quite 
different  from  the  direct  force  at  the  3/4-point,  Theoretically,  the  direct  forces  in  the  arch 
did  not  vary  much  at  these  points,  as  can  be  seen  in  Figure  26.  Close  correlation  between 
the  theoretical  and  experimental  moment  at  the  quarter  points  and  at  the  crown  of  the  arch  is 
evident  in  Figure  27.  Similarly,  the  shearing  force  at  the  location  of  strain  rosette 
SG-13-14-15  is  shown  in  Figure  28.  The  theoretical  and  experimental  load-reaction  curves 
are  compared  in  Figure  29.  Close  correlations  are  evident  in  all  four  reactions. 
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Figure  23.  Comparison  of  theoretical  and  experimental  static  load-deflection  curves 
for  the  antisymmetrically  loaded  arches. 


Arches  AH-2,  AH-3,  and  AH-4  were  preloaded  statically  to  approximately  60  percent  of 
the  yield  load  prior  to  dynamic  testing.  The  data  from  these  test6  are  not  included  in  this 
report;  however,  portions  of  the  data  are  plotted  in  Figures  23  and  29  along  with  the  data  from 
arch  AH-1.  The  permanent  strains  experienced  in  these  static  preload  tests  were  negligible. 

Free-Vibration  Tests 

Efforts  were  made  to  obtain  experimental  values  of  natural  frequencies  of  the  deflection 
mode  and  compression  mode.  However,  only  the  lowest  frequency  of  the  antisymmetrical 
mode  was  obtained.  These  results  are  presented  in  Table  4.  The  average  natural  frequency 
from  the  experimental  data  was  5**  cps,  which  agrees  well  with  the  theoretical  natural 
frequency  of  the  first  antisymmeirical  mode,  60  cps,  as  shown  In  Table  D-5  of  Appendix  D. 
The  theoretical  value  was  obtained  from  a  40-bar  distributed  mass  system. 
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When  viscous  damping  is  present  in  an  oscillating  structural  member,  the  degree  of 
damping  or  the  rate  of  decay  is  measured  by  means  of  the  logarithmic  decrement,  6,  which  is 
expressed  as  follows: 


i 


ikr 


where  is  the  initial  amplitude  and  xn  is  the  amplitude  at  the  n  cycle  from  the  initial  point. 
When  dry  friction  or  Coulomb  damping  is  also  present  with  the  viscous  damping,  the  use  of  the 
above  equation  will  lead  to  erroneous  results.  A  semigraphlcal  trial -and-error  method  has 
been  presented  by  Jacobson  and  Ayre;*“  however,  the  procedure  requires  skilled  Judgment 
and  can  lead  to  inconsistent  results.  A  method  presented  by  DaDeppo11  applies  least-square 
Iteration  to  determine  the  damping  constants  for  a  dry-friction,  viscous-damped  system.  The 
method  is  not  well  suited  to  manual  computation.  A  program  is  presented  in  Appendix  B  for 
calculating  the  damping  constants  by  means  of  a  digital  computer.  The  experimental  decay 
records  of  the  arches  were  analyzed  by  this  computer  program.  The  value  of  the  logarithmic 
decrement,  j,  the  damping  ratio,  ce/cc,  and  the  Coulomb  damping,  A,  are  presented  in 
Table  4.  The  damping  ratio  is  defined  as  the  ratio  of  the  coefficient  of  viscous  damping,  ce,  to 
the  critical  damping  coefficient,  cc.  The  average  value  of  the  damping  ratio  for  the  first  anti- 
symmetrical  mode  was  0.095.  In  the  theoretical  analysis,  the  damping  factor  was  expressed  in 
the  compression  mode,  and  a  corresponding  value  of  0.018  was  used  as  discussed  in  Appendix 
C.  A  typical  free-vlbration  record  is  shown  in  Figure  9. 


Siroin  (  1,000  pin. /in.) 

Figure  25.  Comparison  of  theoretical  and  experimental  strains  at  the  3/4 -point 
of  the  antisymmetrically  loaded  arch. 
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Comparison  of  theoretical  and  experimental  direct  forces  in  the  antisymmetrically  loaded  arch. 
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Comparison  of  theoretical  and  experimental  moment  in  the  antisymmetrically  loaded  arch. 
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theoretical 


Figure  28,  Comparison  of  theoretical  and  experimental  shearing  forces  in  the 
antisymmetrically  loaded  arch  at  strain  rosette  SG-13,  14,  15. 


Dynamic  Tests 

After  a  maximum  static  load  of  243  lb/in.  was  applied  and  unloaded,  arch  AH-2  was 
subjected  to  six  dynamic  loads  of  increasing  magnitudes,  with  peak  loads  between  116  lb/in. 
and  780  lb/in.  The  postshot  view  of  the  arch  is  shown  in  Figure  F-4.  Maximum  values  of 
strains,  reactions,  deflection,  and  acceleration  at  the  3/4 -point  were  obtained  from  data 
recorded  on  the  oscillograph  paper  as  shown  in  Figure  10. 

Similarly,  arches  AH-3  and  AH-4  were  tested  dynamically  to  peak  loads  of  780  lb/ln. 
and  664  lb/in. ,  respectively.  Reduced  data  of  all  the  dynamic  tests  are  presented  in 
Tables  A-4  and  A-5. 

In  comparing  the  strain  values,  it  was  found  that  the  strains  at  the  lower  extremity  of  the 
3/4 -point  were  higher  than  the  same  location  of  the  1/4 -point  for  each  dynamic  load.  At  a 
peak  load  of  780  lb/in. ,  strain  gages  SG-4  and  SG-12  indicated  strains  of  7,090  /uin./in.  and 
6,940  A*tn./ln. ,  respectively;  a  similar  relationship  was  observed  in  the  static  test  data. 

It  should  be  noted  that  certain  dynamic  tests  did  not  yield  valid  data  and  are  therefore  not 
recorded  in  Tables  A-4  and  A-5.  The  strain  data  of  Table  A-4  were  taken  at  each  peak  value, 

which  presented  a  problem  in  correlation  of  theoretical  and  experimental  response  because  all  *  « 

the  peak  strains  did  not  occur  at  the  same  time.  Only  a  portion  of  the  data  in  the  elastic  range 

of  arch  AH-2  was  analyzed,  and  the  resulting  moments,  direct  forces,  principal  and  shearing 

stresses  are  summarized  in  Table  A- 8  and  Figures  A -3  and  A-4.  The  strain  data  of  arches 

AH-2,  AH-3,  and  AH-4  are  plotted  in  Figures  A-5,  A-6,  and  A-7,  respectively. 

The  load-deflection  curves  for  both  theoretical  and  experimental  results  are  given  in 
Figure  30.  The  experimental  load-deflection  curve  represents  the  mean  values  of  the  data 
taken  from  the  maximum  values  of  each  test  of  arches  AH-2,  AH-3,  and  AH-4.  The  theoretical 
curve  is  plotted  from  the  data  of  Appendix  C.  Good  correlation  between  the  theoretical  and 
experimental  curves  is  evident  in  Figure  30. 


31 


-rfi» 


Comparison  of  theoretical  and  experimental  reactions  of  the  antisymmetrically  loaded  arch  in  the  elastic  range. 
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Figure  30.  Comparison  of  dynamic  load-deflection  curves  for 
the  antisymmetrically  loaded  arch. 


Figure  31  3hows  the  comparison  of  the  theoretical  and  experimental  moments  of  the 
crown  and  quarter  points  ot  arch  AH-2.  Only  experimental  moments  in  the  elastic  range  are 
presented  for  comparison.  Moments  at  the  quarter  points  show  good  agreement,  but  the 
moments  at  the  crown  do  not  agree  very  well.  It  is  believed  that  the  disagreement  of  moments 
at  the  crown  of  the  arch  is  due  primarily  to  the  inaccuracy  of  the  small  strain  readings  of 
SG-5  through  SG-8.  The  experimental  values  of  direct  forces  did  not  agree  with  the  theoretical 
values  because  the  peak  values  of  the  strains  did  not  occur  simultaneously.  In  the  calculation 
of  experimental  direct  forces,  the  algebraic  sum  of  different  strains  at  the  upper  and  lower 
extremities  was  taken.  Small  inaccuracies  will  result  in  large  differences  in  the  direct 
force. 

The  axial  forces  obtained  from  the  simplified  analysis  (Appendix  E)  were  very  close  to 
those  obtained  from  the  more  rigorous  analysis  of  Appendix  C.  However,  the  values  of 
bending  moment  were  notably  higher  in  the  same  comparison,  since  the  damping  was  not 
considered  in  the  simplified  analysis.  The  axial  forces  have  time  variations  which  correspond 
to  the  second  symmetrical  mode  and  the  moment  has  a  time  variation  which  corresponds  to  the 
first  antisymmetrical  mode.  Obviously,  the  damping  effect  on  the  moment  was  much  higher 
than  on  the  axial  force. 
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Figure  31.  Comparison  of  theoretical  and  experimental  dynamic  moments. 


GENERAL  DISCUSSION 

Design,  Behavior,  and  Analysis  of  the  Arches 

The  depth  and  width  of  the  arches  were  limited  to  4  inches  and  7-3/4  inches, 
respectively,  because  of  the  physical  dimensions  and  safe-load  capacity  of  the  blast  simulator. 
To  provide  symmetry  of  the  cross  section  with  respect  to  the  neutral  axis  would  have  required 
the  use  of  a  thin-gaged  steel  member.  As  a  compromise,  a  top-heavy  section  was  selected 
for  the  arches,  consisting  of  a  4M13  cold- rolled  section  with  steel  strips,  5/8  inch  by 
1-7/8  inches,  welded  to  each  side  of  the  top  flange.  In  the  cold-roll  process,  the  4M13 
section  was  strained  beyond  yield.  The  flange  strips  were  cold-rolled  separately  and  then 
welded  to  the  4M13  section;  as  a  result  the  flange  strips  were  inelastically  strained  less  than 
the  4M13  section.  The  strain-hardening  properties  of  the  flange  strips  and  the  4M13  section 
were  determined  by  static  and  dynamic  tensile  coupon  tests. 

Because  of  the  strain  hardening  encountered  in  the  cold- roll  process,  the  static  and 
dynamic  stress-strain  curves  of  the  material  did  not  exhibit  a  definite  yield  point.  It  was 
also  observed  that  the  total  elongation  of  the  tensile  specimens  was  less  for  strain-hardened 
steel  than  for  similar  material  that  was  not  strain  hardened.  It  is  believed  that  the  arches 
would  have  deflected  more  at  ultimate  load  if  these  strain-hardening  effects  had  not  been 
present.  Furthermore,  the  ultimate  ioad  of  the  arch  was  limited  by  the  ultimate  strength  of 
the  lower  flange.  As  the  lower  flange  was  stressed  into  the  inelastic  range,  the  neutral  axis 
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moved  closer  to  the  upper  flange  which  caused  the  strain  in  the  lower  flange  to  increase  more 
rapidly.  As  the  arch  approached  the  ultimate  load,  a  considerable  portion  of  the  arch 
remained  in  a  low-stress  state  at  the  critical  section;  thus,  failure  occurred  before  the  full 
ductile  potential  of  the  arch  was  developed. 

In  the  analysis,  it  was  necessary  to  consider  the  nonhomogeneous  section  created  by  the 
4M13  section  and  the  flange  strips.  It  was  desirable  to  express  the  stiffness  properties  of  the 
cross  section  the  products  El  and  EA.  Separation  of  the  modulus,  E,  from  these  product 
terms  wouiq  have  complicated  the  analysis. 

Due  to  the  lack  of  a  definite  yield  plateau  in  the  steel  used  in  the  arches  (particularly  the 
flanges),  a  trilinear  stress-strain  curve  was  used  for  the  static  analysis.  Ordinarily  this 
would  not  be  necessary  since  the  cold-roll  process  is  not  generally  used  in  the  fabrication  of 
large  arches.  The  analytic  method  presented  in  Appendix  B  is  valid  for  materials  with 
bilinear  or  trilinear  stress-strain  curves  or  even  for  a  material  that  does  not  have  an  elastic 
range. 

In  the  elastic  range,  the  static  analysis  initially  involved  solving  a  statically  determinate 
arch  without  the  redundant  reactions  (the  two-hinged  arch  has  one  redundant  reaction  and  a 
fixed-end  arch  has  three).  The  redundant  reaction  can  then  be  obtained  by  satisfying  the 
boundary  conditions  at  the  supports.  In  this  analysis,  elastic  supports  as  well  as  rigid 
supports  can  be  handled. 

In  the  ineiastic  range,  the  incremental  stress  method  was  used  along  with  an  iteration 
process  to  determine  various  locations  of  yielding  until  the  ultimate  strength  was  reached  at 
some  location  in  the  arch.  The  analysis  was  carried  only  as  far  as  the  development  of  the 
ultimate  strength. 

Analysis  of  the  two-hinged  arch  of  uniform  cross  section  indicated  that  moments  at  the 
crown  were  small  for  a  uniform  load  distributed  over  one-half  the  arch.  This  suggested  that 
the  placement  of  a  third  hinge  at  the  crown  would  not  substantially  alter  the  behavior  of  the 
arch.  Thus,  a  three-hinged  arch  would  support  the  loads  as  well  as  the  two-hinged  arch.  The 
analysis  of  a  three-hinged  arch  is  much  simpler  than  that  of  a  two-hinged  arch.  In  addition, 
such  a  statically  determinate  arch  would  not  be  subjected  to  the  effects  of  thermal  stress, 
differential  settling  of  the  support,  and  initial  misalignment.  Furthermore,  the  variation  in 
cross  section  could  be  handled  with  less  difficulty. 

Dynamic  analysis  of  the  arch  employed  a  discrete  framework  in  both  the  elastic  and 
inelastic  ranges.  The  mass  of  the  arch  was  concentrated  at  the  flexible  joints,  which  were 
connected  by  rigid  bars.  An  iteration  process  was  used  to  obtain  the  forces  and  displacements 
at  each  joint  by  satisfying  the  respective  equations  of  motion.  The  numerical  integration 
method  by  Newmark  was  an  essential  tool  in  this  analysis. 12 

An  additional  assumption  was  made  for  the  dynamic  analysis  in  the  inelastic  r  ange.  The 
cross  section  of  the  arch  was  represented  only  by  the  two  flanges.  This  assumption  was 
desirable  due  to  the  complexity  of  the  bilinear  stress-strain  relationship  (Figure  C-5)  in  the 
loading  and  unloading  process.  Accuracy  of  the  dynamic  analysis  by  the  iteration  process 
depends  largely  on  the  simplicity  of  the  idealized  stress-strain  curve  and  the  number  of  the 
joints  in  the  discrete  framework.  Convergence  of  iteration  is  slow  if  these  factors  are 
complicated. 

In  the  determination  of  natural  frequencies  and  mode  shapes  (Appendix  D),  it  was  found 
that  the  distributed-mass  system  provided  better  results  than  the  concentrated-mass  system 
for  the  same  number  of  joints.  However,  it  is  believed  that  the  use  of  a  distributed-mass 
system  for  the  dynamic  response  theory  would  undesirably  complicate  the  analysis  and 
increase  calculation  time. 

In  the  actual  arch  tests,  there  were  numerous  contributions  to  total  damping  of  the 
system:  viscous  damping  (internal  and  external).  Coulomb  damping,  and  damping  at  the 
supports.  To  separate  the  effect  of  each  type  of  damping  would  require  extensive  study  of  the 
properties  of  the  material;  this  study  was  not  justified  from  an  engineering  standpoint.  The 
damping  factor  used  in  the  theoretical  "evelopmenl  was  obtained  experimentally  from  the 
free-vibration  tests  in  the  deflection  mode  and  represented  the  combined  effect  from  all 
sources.  It  is  believed  that  the  supports  provided  a  major  contribution  to  the  total  damping; 
the  number  and  quality  of  the  support  connections  greatly  influenced  the  damping  charac¬ 
teristics  of  the  arch.  Also,  it  is  believed  that  a  three-hinged  arch  would  have  higher  damping 
than  a  two-hinged  arch  or  a  hingeless  arch,  provided  other  factors  remained  unchanged. 
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The  simplified  dynamic  analysis  of  Appendix  E  is  a  rapid,  easy  method  of  obtaining 
information  on  a  portion  of  the  dynamic  response  of  an  antisymmetrically  loaded  arch.  This 
method  is  limited  to  the  response  erf  the  quarter  point  of  the  loaded  side  of  the  arch.  The 
damping  effect  was  not  considered  in  the  analysis  and,  as  a  result,  the  moment  values  were 
larger  than  their  counterparts  in  Appendix  C.  The  values  of  axial  forces  obtained  from  the 
simplified  analysis  compared  very  well  with  those  of  the  more  rigorous  analysis. 

The  digital  computer  proved  indispensable  in  both  the  static  and  dynamic  analyses  of  the 
arch  system,  particularly  when  elastic  supports  were  involved. 

Effect  of  Load  History 

The  effect  of  load  history  on  the  dynamic  response  of  structures  had  been  reported  for 
structural  connections  and  portal  frames.  1*^|3  in  the  static  test  of  arch  AH-1,  pressure 
leakage  necessitated  unloading  and  reloading  the  arch  after  certain  portions  of  the  lower  flange 
had  been  strain-hardened.  The  load -deflection  curve  for  the  reloading  was  higher,  before  it 
became  flat,  than  the  initial  loading  curve.  This  effect  is  known  as  strain  aging.  5  When  the 
metal  is  strained  and  then  aged,  the  sharp  yield  and  the  yield  elongation  may  return  and 
increase  the  resistance  of  the  material.  This  phenomenon  is  accompanied  by  a  decrease  in  the 
ductility  of  the  steel.  5  Multiple  loading  of  the  arch  will  reduce  the  flexibility  of  the  arch  in 
the  inelastic  range,  particularly  near  the  ultimate  load.  The  cold-roll  fabrication  process  of 
the  arch  can  be  considered  a  form  of  loading;  its  effects  have  been  discussed  in  the  previous 
section. 

Application  of  Results 

Correlation  of  the  theoretical  and  experimental  data  provides  a  basis  for  design  of  arch 
structures  to  resist  blast  loading.  The  information  contained  in  this  report  can  be  used  to 
refine  the  design  and  fabrication  of  arch  structures. 

This  study  applies  only  to  the  loading,  support  conditions,  and  configuration  of  the  arch. 
Generalization  of  the  results  to  produce  a  simplified  method  for  the  solution  of  arch  problems 
should  be  performed  with  caution,  hi  relating  the  results  of  this  investigation  to  design 
practice,  particular  care  must  be  taken  to  examine  the  load  and  mass  distributions,  stiffness 
characteristics,  and  support  condition  of  the  structure.  The  methods  and  computer  program 
for  the  theoretical  solution  of  arch  behavior  provided  in  this  report  can  be  used  to  study  arches 
of  other  configurations  and  loading,  provided  proper  input  parameters  are  employed. 

Accuracy  of  Measurements 

The  instrumentation  used  in  this  experiment  was  carefully  selected  and  calibrated  before 
testing.  The  obtained  measurements  are  believed  to  be  accurate  within  the  manufacturer- 
recommended  tolerance  of  each  measuring  device.  Data  reduction,  although  subject  to  human 
error,  was  carefully  performed  and  checked  to  eliminate  all  significant  errors.  The  reduced 
experimental  data  are  assumed  to  have  the  following  accuracy: 

Data  Accuracy  (%) 


Static: 


Pressure 

2  to 

3 

Reaction 

2  to 

3 

Deflection 

2  to 

3 

Strain 

4  to 

5 

Dynamic: 

Pressure 

4  to 

5 

Reaction 

4  to 

5 

Deflection 

4  to 

5 

Strain 

8  to 

10 

Strain  rate 

10  to 

12 

Acceleration 

12  to 

15 
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Conversion  of  the  experimental  data  for  comparison  with  the  theoretical  solutions  was 
performed  by  a  computer.  The  results  are  believed  to  be  as  accurate  as  the  assumptions 
permit. 

FINDINGS  AND  CONCLUSIONS 

The  experimental  tests  and  the  theoretical  solutions  included  in  this  report  are  limited  to 
a  specific  two-hinged  arch  with  certain  loading  conditions;  however,  some  of  the  findings  and 
conclusions  are  applicable  to  arches  with  other  loading  and  support  conditions  and 
configurations. 

General 

1.  Tensile  coupon  tests  Indicated  that  the  steel  used  in  the  arches  was  strain  hardened;  the 
material  lacked  a  definite  yield  plateau  common  to  structural  steel  and  consequently  was  better 
represented  by  a  trilinear  stress-strain  curve  than  by  the  usual  bilinear  curve. 

2.  The  ductility  of  the  arch  material  was  reduced  primarily  by  the  strain  resulting  from  the 
cold-roll  process;  the  arch  failed  before  the  ductile  potential  was  developed. 

3.  When  the  material  was  strain  hardened  and  then  aged,  as  when  the  arch  was  unloaded  and 
loaded  in  the  Inelastic  range,  the  stress  rose  above  the  original  stress-strain  curve  and  the 
ductility  of  the  material  was  considerably  reduced. 

4.  The  failure  modes  for  both  the  statically  and  dynamically  loaded  arches  were  similar.  The 
governing  factor  was  the  ultimate  stress  at  the  lower  flanges  of  the  quarter  points. 

Static  Tests 

1.  Good  correlation  was  obtained  between  experimental  behavior  and  the  theoretical  analysis 
in  the  elastic  and  inelastic  ranges. 

2.  The  deflected  shape  of  the  arch  was  basically  antisymmetries!,  with  a  rotation  of  the  chord 
connecting  the  two  hinges  introduced  by  the  yielding  supports. 

3.  Bending  moment  was  the  predominant  factor  in  the  stress  and  strain  distribution  in  the 
arch;  the  direct  axial  forces  exerted  only  a  minor  effect  in  developing  the  stress  and  strain  in 
the  arch. 

4.  The  yielding  support  on  the  loaded  end  reduced  the  axial  forces  in  the  arch,  increased  the 
bending  moment  and  deflection  on  the  loaded  half  of  the  arch,  and  decreased  the  moment  and 
deflection  on  the  unloaded  half  of  the  arch. 

5.  The  maximum  bending  moment  and  deflection  occurred  very  near  the  quarter  points  of  the 
arch,  where  the  shearing  forces  were  zero. 

6.  The  maximum  stresses  occurred  in  the  lower  flanges  of  the  quarter  points. 

7.  The  crown  of  the  arch,  which  was  very  close  to  the  point  of  contraflexure,  maintained  a 
very  low  bending  moment  throughout  the  loading  history. 

8.  Collapse  occurred  when  the  ultimate  stress  at  the  lower  flanges  of  the  quarter  points  was 
reached. 

Dynamic  Tests 

1,  The  theoretical  analyses  satisfactorily  predicted  the  experimental  dynamic  behavior  of  the 
arches  in  the  elastic  and  inelastic  ranges.  The  simplified  analysis  proved  to  be  sufficiently 
accurate  for  design  purposes. 


2.  The  time-dependent  displacements  occurred  at  a  frequency  corresponding  to  the  first 
antisymmetrical  mode  of  the  free  vibration  of  the  arch.  This  indicated  that  the  first 
antisymmet.  ical  mode  was  predominant  in  the  response. 

3.  The  axial  force  remained  essentially  uniform  along  the  length  of  the  arch  and  varied  with 
time  according  to  the  frequency  of  the  second  symmetrical  mode  of  the  arch  (breathing  mode). 

4.  The  bending  moments  varied  with  time  at  the  frequency  of  the  first  antisymmetrical  mode 
of  the  arch.  Although  some  contribution  of  the  third  antisymmetrical  mode  was  observed,  the 
first  antisymmetrical  mode  definitely  made  the  main  contribution  to  the  moment. 

5.  Maximum  stresses  occurx'ed  at  the  lower  flanges  of  the  quarter  points  primarily  due  to 
bending  action. 

6.  Bending  moment  at  the  crown  of  the  arch  was  very  small  for  all  dynamic  tests. 

7.  Damping  of  the  system  had  an  important  effect  on  the  bending  response  of  the  arch;  its 
effect  on  the  axial  force  was  small. 

DESIGN  RECOMMENDATIONS 

The  methods  used  in  the  design  of  arch  structures  depend  on  such  factors  as  the  material 
used,  geometry  of  the  arch,  conditions  of  loading  and  support;  and  the  distribution  of  mass 
and  stiffness  along  the  arch.  The  two-hinged  arch  discussed  in  this  report  is  only  one  specific 
case.  However,  certain  design  recommendations  can  be  drawn  from  the  theoretical  analyses 
and  experimental  data  which  apply  to  two-hinged  circular  arches  subjected  to  overpressure 
from  a  traveling  blast  wave. 

Static  Design  and  Analysis 

The  ultimate  strength  of  the  antisymmetrically  loaded  two-hinged  arch  is  essentially 
limited  by  the  moment-resisting  capacities  at  the  quarter  points  of  the  arch.  This  suggests 
that  special  reinforcement  of  the  arch  sections  in  the  form  of  cover  plates,  or  increased  depth 
at  the  critical  sections  will  considerably  increase  the  ultimate  strength  at  a  small  increase  in 
fabrication  cost. 

The  bending  moment  at  the  crown  remained  low  throughout  the  entire  loading  history. 
This  suggests  tuat  a  three-hinged  arch  (hinges  at  the  supports  and  die  crown)  of  the  same 
cross  section  will  respond  very  much  like  the  two-hinged  arch.  Thus  it  would  be  logical  to 
consider  a  three-hinged  arch  as  the  first  alternate  design  of  the  two-hinged  arch.  As  a 
statically  determinate  stiucture,  the  three-hinged  arch  is  less  affected  by  thermal  stress, 
differential  settling  of  supports,  and  initial  misalignment.  The  stress  level  is  more  readily 
controlled  for  certain  loading  conditions. 

The  arch  response  is  very  sensitive  to  the  support  conditions.  It  is  essential  that  proper 
eva.uatlon  of  the  structural  foundation  be  made  before  designing  the  arch.  For  arches  resting 
on  elastic  supports,  the  axial  forces  in  the  arch  will  decrease  and  the  bending  moments  will 
increase  with  increasing  flexibility  of  the  supports.  In  case  of  uncertainty,  the  design  should 
be  adequate  for  the  twe  worst  conditions. 

The  imbalanced  arch  cross  section  has  a  tendency  to  limit  the  moment-resisting  capacity 
of  the  arch  according  to  the  yield  strength  of  one  flange.  As  the  inelastic  strain  develops,  the 
neutral  axis  moves  away  from  the  flange  having  higher  stress  and  thus  accelerates  the 
straining  process.  For  isotropic  materials  (most  metals  can  be  considered  isotropic  for 
design  purposes),  a  balanced  cross  section  is  essential  for  the  most  economical  design. 

Strain  hardening  of  the  steel  reduces  its  ductility;  as  a  result,  collapse  of  a 
strain -hardened  arch  occurs  at  the  ultimate  load  with  little  warning.  The  use  of  strain- 
hardened  material  is  not  recommended  for  arch  construct  .on.  In  particular,  a  cold-roll 
process  should  not  be  used  to  fabricate  arches  other  than  model  arches  and  sheet  metal 
arches. 


38 


The  use  of  a  digital  computer  can  be  of  great  value  In  the  arch  analysis  particularly  when 
elastic  supports  are  involved.  The  programs  presented  in  Appendix  B  may  be  readily  modified 
for  arches  of  different  geometries  and  support  conditions. 

Qiasi-pseudo  static  response  may  be  obtained  from  the  dynamic  solution  by  introducing  a 
high  level  of  damping.  A  damping  factor  of  0, 3  proved  to  be  satisfactory  for  the  arch  under 
consideration. 

Dynamic  Design  and  Analysis 

For  long-duration  loads  (in  which  the  ratio  of  the  load  duration  to  the  natural  period  of 
the  arch  is  greater  than  ten),  dynamic  bending  moments  at  the  quarter  points  are  equal  in 
magnitude.  However,  the  combination  of  moment  and  direct  force  will  cause  the  lower  flange 
of  the  quarter  point  on  the  loaded  portion  of  the  arch  to  yield  first. 

Maximum  stresses  are  contributed  primarily  by  bending.  Special  reinforcement  by 
varying  the  effective  depth  of  the  arch  sections  at  the  quarter  points  is  desirable.  For 
dynamically  loaded  structures,  the  increased  depth  is  more  efficient  reinforcement  than  cover 
plates  because  less  mass  is  involved. 

Rapidly  strained  steel  exhibits  a  considerable  increase  in  dynamic  yield  point;  references 
2  and  3  report  results  for  initially  unstrained  specimens  of  A-7  mild  steel.  No  research  has 
been  found  which  reports  the  increase  in  dynamic  strength  of  steel  after  it  had  been  previously 
strained  into  the  strain- hardening  region.  Indications  from  this  investigation  are  that  the 
percentage  increase  in  dynamic  strength  properties  diminishes  after  previous  inelastic 
straining  has  occurred.  Specimens  from  some  of  these  arches  exhibit  strength  increases  of 
about  10  percent  for  strain  rates  of  0.5  in./in./sec. 

The  level  of  damping  is  very  important  for  the  calculation  of  dynamic  response  of  the 
arch.  A  damping  factor  of  0. 1  will  probably  be  appropriate  for  most  steel  arches. 

Due  to  the  absence  of  the  damping  terms,  the  use  of  the  simplified  analysis  in  Appendix 
E  is  limited  to  systems  of  light  damping.  This  is  demonstrated  in  the  comparisons  of  axial 
force  and  moment  obtained  from  the  simplified  and  more  rigorous  dynamic  solutions. 

Similar  to  the  static  response,  the  bending  moment  at  the  crown  for  dynamic  load 
remains  low.  Consideration  of  a  three-hinged  arch  is  also  valid  for  dynamically  loaded 
arches.  In  addition,  the  damping  value  will  be  higher  for  a  three-hinged  arch  of  similar 
cross- section  properties. 

The  efficient  use  of  a  digital  computer  is  essential  in  obtaining  accurate  dynamic 
response  of  the  arch.  The  dynamic  solution  by  the  method  of  numerical  integration- is 
impractical,  if  not  impossible,  by  manual  calculation. 
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LIST  OF  SYMBOLS 

General  Note: 

Superscripts:  Refer  to: 

f,  r  left  and  right  sides  of  a  joint 

x  first  derivative  of  x  with  respect  to  time 

x  second  derivative  of  x  with  respect  to  time 
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Subscripts: 
h,  i,  j,  k,  i,  m 

x.  y 
r,  t 
u,  * 
t,  b 


designations  for  location  of  joints 
designations  for  x,  y  direction 
radial  and  tangential  components 
upper  and  lower  flanges 

top  and  bottom  flanges  of  the  two-flange  representation  of  the 
cross  section 


A  symmetrical  matrix  for  natural  mode  and  buckling  load  equation 

A^m  static  stress  amplification  factor  at  joint  j 

Am  static  amplification  factor  for  the  simplified  analysis 

A.,  A.  area  of  top  and  bottom  flanges  of  the  two-flange  representation  of  the 
cross  section  (in. ) 

B  symmetrical  matrix  for  natural  mode  and  buckling  load  equation 

B'  symmetrical  matrix  for  natural  mode  and  buckling  load  equation 

b  instantaneous  set 

C  constant 

C  symmetrical  matrix  for  natural  mode  and  buckling  load  equation 

Cm  initial  misalignment  of  arch  supports  (in. ) 

c  parameter  used  in  the  equations  for  natural  modes  and  buckling  loads 

cc  critical  damping  coefficient 

ce  external  damping  coefficient 

c  ,  Cj  distances  from  neutral  axis  to  upper  and  lower  flanges  (in. ) 

D  (x)  dynamic  amplification  factor  for  simplified  analysis 

d  depth  of  arch  cross  section 

E  modulus  of  elasticity  (psi) 

E'  inelastic  modulus  of  material  (psi) 

Ej  second  modulus  of  material  (psi) 

Eg  third  modulus  of  material  (psi) 

EAj  extensional  stiffness  parameter  of  joint  j  (ib) 

o 

Elj  flexural  stiffness  parameter  of  joint  j  (lb-in.  ) 
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e 

F 


f 

fn 

fy 

H 

k 

k 

L 


M, 


M 


m, 


N, 


n 

n 


pjx’  pJy 

*V  PJy 

PJr*  Pit 
P 
R 


“1 

Tn 

T(x) 


W 


E 

WT 


Wjr’  WJt 


w 


strain  (in.  /in. ) 
external  force  (lb) 
stress  (psi) 

natural  frequency  of  vibration  (cps) 

yield  stress  (psi) 

horizontal  redundant  reaction  (lb) 

parameter  used  in  the  equations  for  natural  modes  and  buckling  loads 

spring  constant  of  elastic  support  (in./lb) 

length  of  bars  for  the  discrete  framework  (in. ) 

span  of  the  arch  (in. ) 

bending  moment  at  joint  j  (in.  -lb) 

plastic  moment  (in.  -lb) 

mass  of  joint  j  (lb-secfyin.  2) 

axial  force  due  to  the  deformation  of  bar  j  (lb) 

number  of  bars  for  the  discrete  framework  (in. ) 

a  constant  used  in  the  simplified  analysis 

x,  y  components  of  external  forces  acting  at  joint  j  (lb) 

x,  y  components  of  internal  forces  at  joint  j  (lb) 

radial  and  tangential  forces  at  joint  j  (lb) 

pressure  (psi) 

radius  of  the  arch  (in. ) 

reaction  1,  etc.  (lb) 

natural  period  (msec/cycle) 

natural  periods  used  in  the  simplified  analysis  (msec/cycle) 
time  (sec) 

duration  of  dynamic  load  (sec) 
total  external  work  (in.  -lb) 
total  internal  work  (in.  -lb) 

external  resultant  forces  at  joint  j  in  radial  and  tangential  directions  (lb) 
uniform  load  (Ib/in. ) 
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w  buckling  load  used  in  the  simplified  analysis  (lb/in. ) 

Z  plastic  section  modulus  (in.  3) 

Zjr,  Zjt  internal  resultant  forces  at  joint  j  in  radial  and  tangential  directions  (lb) 
a  angle  between  radii  of  support  and  a  point  on  the  arch  (rad) 
ft  angle  between  radial  of  a  support  and  x  axis  (rad) 
coefficient  of  numerical  integration 
y  maximum  shearing  stress  in  the  web 
6  logarithmic  decrement 

dj  deformation  of  bar  j  (in. ) 

dj  first-order  approximation  of  d.  (in. ) 

A  deflection  (in. ) 

djx>  djy  x,  y  components  of  displacement  at  joint  j  (in. ) 
djr,  d<t  radial  and  tangential  displacements  at  joint  j  (in. ) 
dt  time  Increment  (sec) 

dj  angular  change,  elastic  weight  at  joint  j  (rad) 

first  order  approximation  of  8^  (rad) 

0. '  inelastic  angular  change  at  joint  j  (rad) 

X  parameter  used  in  the  equations  for  natural  modes  and  buckling  loads 

2  3 

H  mass  per  unit  length  of  the  arch  (lb-sec  /in.  ) 
v  damping  factor 

v'  damping  factor  associated  with  "breathing"  mode  of  vibration 

p  parameter  used  in  the  equations  for  natural  modes  and  buckling  loads 
a j  ^  principal  stresses  1  and  2 

4>q  central  angle  of  the  arch  (rad  or  deg) 

<t> ^  central  angle  subtended  by  bar  j  (rad) 

<£p  angle  of  principal  axes  with  respect  to  strain  rosette 
ip  angular  rotation  (rad) 

ip  j  angular  rotation  of  bar  j  (rad) 

ip  I  first  order  approximation  of  ip^  (rad) 
u  circular  frequency  of  vibration  (rad/sec) 
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Appendix  A 

MEASUREMENTS  OF  PRESSURE,  DEFLECTION,  AND  STRAIN 


Table  A-l.  Static  Test  Data  for  Antisymmetrically  Loaded  Arch  AH-1,  Test  No.  2 


Table  A-3.  Static  Test  Data  for  Antisymmetrically  Loaded  Arch  AH-1,  Test  No 
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Table  A-6.  Summary  of  Maximum  Dynamic  Data  for  Arch  AH-2 


Load 
(lb/in. ) 

SG-1 

(//in.  /in. ) 

SG-4 

(//in.  /in. ) 

Moment 
(in.  -kips) 

Direct 

Force 

(kips) 

116 

-460 

810 

129.0 

3. 29 

235 

-800 

1,430 

227.0 

14.10 

411 

-1,480 

2,510 

406.  0 

34.50 

SG-5 

SG-8 

116 

-170 

-380 

-20.8 

42.60 

235 

-270 

-460 

-18.5 

59.  00 

411 

-470 

-800 

-33.8 

SG-9 

SG-12 

116 

370 

-760 

235 

590 

-1,260 

411 

1, 140 

-2, 200 

v-  1 

-10. 20 

Load 
(lb/in. ) 

Principal  Stress 

Phi 

(rad) 

Max  Shear 
Stress 
(ksi) 

Rad  Shear 
(kips) 

(ksi) 

a 2  (ksi) 

116 

4.63 

7.44 

-6. 15 

235 

8.64 

12.90 

411 

17.  20 

■ 

-0.255 

21.60 

-19.40 

i  ■ 

i 
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■ 


m 


Load  ( Ib/in. 


Direct  Force,  upper  scale  (kips).  Moment,  iotver  scale  (in. -kips) 


Figure  A-3.  Maximum  values  of  experimental  moment  and  direct  force  for 
dynamically  leaded  arch  AH-2. 


Figure  A-4.  Maximum  values  of  experimental  shear  and  principal  stress  for  dynamically  loaded  arch  AH-2  at  SG-13,  14,  15. 


SG-16  SG-14 

SG-15  \  I  SG- 


Figure  A-5.  Maximum  dynamic  strains  for  antisymmetrically  loaded  arch  AH-2. 


Maximum  dynamic  strains  for  antisymmetrically  loaded  arch  AH-3. 


CO 


Maximum  dynamic  strains  for  antisymmetrically  loaded  arch  AH-4. 
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Appendix  B 

STATIC  ANALYSIS  OF  TWO-HINGED  CIRCULAR  ARCH 


General 

The  static  analysis  of  the  two-hinged  circular  arch  with  a  uniform  load,  w,  distributed 
over  one-half  the  span  (Figure  B-l),  includes  theoretical  solutions  in  both  the  elastic  and 
inelastic  ranges.  To  simplify  the  numerical  solution,  an  idealized  arch  (Figure  B-2)  is 
taken  along  the  centerline  of  the  actual  arch.  The  elastic  characteristics  are  represented 
by  the  flexible  joints  connected  by  rigid  bars.  ^  This  substitute  framework  is  shown  in 
Figure  B-3. 

Elastic  Range 

In  determining  the  theoretical  behavior  of  the  circular  arch  in  the  elastic  range,  the 
following  assumptions  were  made: 

a.  A  linear  distribution  of  stress  and  strain  is  maintained  at  any  cross  section  of 
the  entire  arch. 

b.  The  arch  cross  sections  are  uniform  between  any  two  adjacent  joints. 

c.  Segments  of  the  arch  are  assumed  to  be  straight. 

d.  Continuity  of  the  arch  is  maintained  in  the  elastic  range. 


Figure  B-l.  Geometry  and  loading  condition  of  the  circular  arch. 
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Figure  B-2.  Coordinate  system  and  loading  condition  for  static  analysis  of  the  circular  arch. 


e.  The  applied  loads  are  assumed  uniform  between  any  two  adjacent  joints  and  can  be 
represented  by  components  of  concentrated  loads. 

f .  Rotations  due  to  elastic  strain  take  place  at  the  joint. 

g.  The  elastic  support  is  represented  by  a  linear  spring. 

h.  The  effect  of  dead  load  is  negligible. 

Assumption  g  is  due  to  the  horizontal  movement  of  the  right  support  during  the  static 
tests.  The  horizontal  permanent  displacements  recorded  were  1  16  inch.  It  is  believed 
that  the  elastic  displacement  during  the  test  was  larger.  The  elastic  constant  was  taken  to 
be  6. 25  x  10  in.  lb  in  the  static  analysis. 

The  two-hinged  arch  considered  is  a  statically  indeterminate  structure  with  a 
redundant  reaction.  If  the  horizontal  restraint  at  n  is  removed,  the  structure  becomes 
readily  solvable  by  the  principles  of  statics.  The  deformed  structure  is  shown  in  Figure  B-4. 
The  simulated  structure,  obtained  by  the  addition  of  the  redundant  reaction.  H,  is  shown  in 
Figure  B-5.  By  applying  the  principle  of  superposition,  the  magnitude  of  H  is  obtained  for 
the  nonyielding  support  by  equating  the  horizontal  deflection  at  n  (dx”  in  Figure  B-5)  to 
the  deflection  Ax'  in  Figure  B-4.  For  structures  with  yielding  supports,  special  boundary 
conditions  must  be  satisfied. 

The  slope -del lection  relationship  of  a  segment  of  the  arch  may  be  expressed  by  the 
moment-area  method. 
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Figure  B-5.  Original  arch  simulated  by  addition  of  redundant  reaction  H. 


For  the  entire  arch,  the  origin  of  the  coordinate  system  is  taken  at  the  left  support 
of  the  arch  as  shown  in  Figure  B-2.  Note  that  the  x  and  y  axes  coincide  with  the  direction 
of  reactions,  R3  and  R4,  respectively.  By  trigonometry,  the  coordina.es  of  any  joint  j  on 
the  arch  may  be  established  as  follows: 


=  R  J^cos  p  (1  -  cos  a)  -  sir.  p  sin  aj 


(B-l) 


y.  =  R  |sin  P  *  sin  (a  -  P )] 


(B-2) 


where  p  is  the  initial  angle  between  the  radial  line  of  the  arch  at  point  0  and  the  horizontal 
axis.  Because  of  symmetry  at  the  crown  of  the  arcii 


(B-3) 


where  <t>Q  is  the  central  angle  of  the  arch. 
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Consider  a  discrete  system  with  a  finite  number  of  joints  connected  by  straight  bars. 

Let  the  angle  between  radii  to  two  adjacent  joints  be  <t>  and  the  central  angle  of  the  arch  be 
</>0\  then 

1 

fc 

.  _  o 
*  -  — 

where  n  is  the  number  of  bars  used  in  the  discrete  system.  The  bar  length  is 

d> 

L  =  2  R  sin  y 

(B-4) 

The  x  and  y  components  of  the  bar  connecting  joints  j  and  i  are 

i- 

X.  =  X.  -  X. 

JX  ]  i 

(B-5) 

>T 

i 

>r 

ii 

>» 

■«— * 

IB-6) 

where  joint  i  is  the  joint  j-1. 

Proceeding  on  the  assumption  that  the  applied  loads  are  uniformly  distributed  along  each 
segment  of  the  arch,  it  is  convenient  to  resolve  the  distributed  loads  into  x  and  y  components 
of  a  concentrated  load  acting  at  the  centroid  of  the  segment.  If  the  segment  is  sufficiently 
small  to  warrant  the  assumption  of  the  member  being  straight  for  the  segment,  the  uniform 
load,  w,  may  be  resolved  into  components  as  follows: 

p.  =  w  v- 

(B-7) 

i 

! 

and 

1 

pjy  =  wxix 

(B-8) 

| 

1 

The  conditions  of  equilibrium  require  that 

i 

i 

i 

i 

n 

v  +  v  +7  wy.  =  u 
nx  ox  jy 

j  =  o 

(B-9) 

i 

i 

n 

i 

I 

i 

p  +  p  -*  y'  wx.  =  o 
ny  oy  jx 

j=o 

(B-10) 

i 

1 

i 
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and 


M 


n 


M  +  i*  y  + 
o  ox  7n 


M 


P  y  + 
ox  Jn 


Poy  %  *  E  [pjx  (y„  -  'j  *  %)  *  p),  fn  -  *i  *  3r)] 
Poy  *n  '  E  *,[1’iy(yn  -  yJ  '  *  xi«(\  -  *)  *  3?)] 

j=0 


(B-ll) 


where  the  subscripts  o  and  n  denote  beginning  and  ending  points  of  the  free  body  considered 
for  equilibrium.  It  in  important  to  note  that  point  j  can  be  an  ending  point  as  well.  A  review 
of  the  statically  determinate  structure  of  Figure  B-4  Indicates  that  the  x  component  of 
reaction  n  is  zero. 


or 


P  =0 
nx 


11 

pox  '  E  wyiy  *  0 


1=0 


n 

P  =  -  wy. 
ox  t—t  'jy 


]=o 


(B-12) 


(B-13) 


The  summation  of  moments  at  joint.  0  is  zero. 


n 

P  .  x  +  /  w  [y.  y.  +  x.  x.  J  = 

ny  n  L-i  V  1  Jy  J  J*/ 


J  =  o 


n 

Pny  s  -f  E  W(yJ  y)y  +  Xj  xix) 

n  j  =o 


(B-14) 


Summing  up  all  forces  in  the  y  direction,  we  obtain 

n 

Poy  ”  "Pny  "  wxjx 
j  =  o 


(B-15) 


Substituting  Equations  B-13  and  B-15  into  Equation  B-ll,  the  bending  moment  at  any  joint 
can  be  calculated. 
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Consider  an  arbitrary  moment  distribution  in  the  vicinity  of  joint  j  as  shown  in 
Figure  B-6,  and  subdivide  the  moment  diagram  into  a  number  of  triangular  prions.  The 
segment  size  iB  chosen  so  that  the  variation  of  moment  can  be  adequately  approximated  by 
straight  lines  between  adjacent  joints.  The  bending  moments  tributary  to  joint  j  can  be 
represented  by  the  triangle  abc  and  the  area  of  the  triangle  is 


6 


j 


ML  L 


El. 


(B-16) 


where  0j  denotes  angular  change  at  j  due  to  elastic  strain  between  dg.  The  areas  abc  and 
defg  are  equivalent,  and  Mj/EI  can  be  taken  as  the  average  elastic  weight  between  d  and  g. 

For  small  angular  changes,  the  relative  displacement  of  joint  j  with  respect  to  any 
point  A,  shown  in  Figure  B-7,  may  be  expressed  in  component  terms  as  follows:  14 


M.  L 


"i  xi  =  “rip 


M.  L 

=  Vi!  Trip 
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Figure  B-7.  Joint  rotation  due  to  eias  ti^  strains. 


Summing  from  0  to  n  and  adding  the  initial  displacement  terms,  the  x-  and 
y-components  of  final  displacement  become 


«  WLL 


A  =  A 
ny  oy 


♦  *  E  -Ari 

j-o  1 


A  =  A 
nx  ox 


JL.  MiL 

+  Vn +  L  ikrh 

i  =  o  1 


<B  17) 


(B-18) 


It  should  be  noted  that  Mj  may  be  obtained  from  Equation  B-ll  by  substituting  j  for  n. 

To  maintain  continuity,  the  sum  of  total  angular  changes  between  0  and  n  must  be 
equal  to  zero  and  Equations  B-17  and  B-18  must  be  satisfied.  This  results  in 


9 . 


j=o  i 


(B-19) 
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For  the  structure  shown  in  Figure  B-4,  dny  *  doy  =  It  follows  from  Equation  B-17 
that 


e 

o 


M.  L 

Tslp 


(B-20) 


Substituting  Equation  B-20  into  B-13  and  noting  that  AQX  =  0 


nx 


M.  L 


n  j  =  0 


.  MiL 

EIJXl  +  S  EIjy 


J 


By  symmetry 


x 


n 


1 


The  above  equation  is  then  reduced  to 


(B-21) 


This  equation  gives  the  horizontal  displacement  of  the  statically  determinate  arch  shown 
in  Figure  B-4. 

In  order  for  this  arch  to  act  as  a  two-hinged  arch,  the  horizontal  restraining  force,  H, 
must  be  provided  as  shown  in  Figure  B-5.  For  simplicity  of  presentation,  let  H  be  unity. 
The  equilibrium  in  the  x  direction  requires  that 


ox 


’  -Pnx  -  1 


(B-22) 


Talcing  moments  about  point  0  and  noting  the  geometric  symmetry  and  equilibrium  condition 
in  the  y  direction,  it  is  clear  that 


i 


(B-23) 


Substituting  Equations  B-22  and  B-23  into  Equation  B-ll  for  new  bending  moment  terms 


M. 

J 


yj 


x. 

J 


(B-24) 
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Substituting  Equation  B-24  into  Equation  B-21  for  horizontal  displacement  in  the  direction  of 
the  unit  load  at  n.  and  introducing  the  redundant  reaction,  H 


A  ' 

x 


(B-25) 


In  the  case  of  a  nonyielding  support 


V  +  V  =  0  (B_26) 


By  substituting  Equations  B-21  and  B-25  into  Equation  B-26,  one  obtains  the  following 
expression  for  the  redundant,  H 


H 


P3 


(B-27) 


For  elastic  support  represented  by  a  linear  spring  with  constant  k  in  the  direction  of  H 


V  +  V  -  kH 

therefore 


(B-28) 


(B-29) 


Note  that  Equation  B-27  is  a  special  case  of  Equation  B-29  with  k  =  0. 

Once  the  redundant,  H,  is  known,  the  two- hinged  arch  can  be  solved  as  a  statically 
determinate  structure  as  outlined  above. 

The  bending  moment  and  the  x-  and  y-components  of  forces  at  each  joint  can  be 
calculated  by  Equations  B-9,  B-10,  and  B-ll.  It  is  desirable  to  convert  the  x-  and 
y-components  of  the  forces  at  each  joint  into  radial  and  tangential  forces  as  follows: 


Fjr  E  pjxeos  {a  ~  fi) 
Pjt  =  Pjy  cos  (a  ' 


‘  jy 


i* 


sin 


sin 


(«  -fi) 

(B-30) 

(a  -0) 

(B-31) 
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where  Pjr  and  Pu  are  the  radial  and  tangential  forces  at  the  Joint  J,  generally  known  as 
radial  shear  and  airect  forces  in  the  arch. 

The  stress  and  strain  of  upper  and  lower  fibers  at  each  section  can  be  calculated  from 
the  bending  moment  and  direct  forces  of  Equations  B-ll  and  B-13.  The  expressions  for 
strains,  eu  and  ej,  for  the  upper  fiber  and  the  lower  fiber  are  .'s  follows: 


e“  ■  ft  *  ft) 


(B-32) 


n  - 


(B-33) 


where  cu  and  C|  are  distances  from  the  neutral  axis  to  the  upper  and  lower  fibers, 
respectively.  The  stresses  fu  and  f j  for  the  upper  and  the  lower  fibers  may  be  expressed 
in  terms  of  strains  as  follows: 


e  E 
u  u 


(B-34) 


(B-35) 


where  Eu  and  Ej  are  the  elastic  moduli  of  the  upper  and  lower  fibers,  respectively. 

The  x-  and  y-componente  of  deflection  can  be  obtained  from  Equations  B-21  and  B-22 
and  can  be  converted  into  radial  and  tangential  components,  djr  and  4jt,  as  follows: 


djr  =  djx  cos  (a  -  0)  +  djy  sin  (a  -  0) 


(B-36) 


^jt  =  4jyco8^a  ~ ’  d.xsbi(«  -  0) 


(B-37) 


Secondary  bending  moments  are  produced  by  direct  forces  acting  on  the  radial  displacements, 
djt,  and  initial  misalignment  of  the  arch.  The  secondary  stresses  are  ordinarily  neglected 
in  structural  analysis;  however,  in  this  arch  analysis  they  are  quite  significant. 


Mr  -  ri.  (■ 


*  cm) 


(B-38) 


Mj'  denotes  the  secondary  bending,  and  Cm  is  the  initial  misalignment.  It  is  evident  that 
the  deflections,  stresses,  and  strains  are  amplified  by  Mj'  as  follows: 


M.  +  M. 


(3-39) 


4 


The  amplification  factor,  Ajm,  should  be  used  as  a  coefficient  in  Equations  B-36,  B-37, 
B-ll,  B-32,  B-33,  B-34,  and  B-35  for  the  final  deflections,  moment,  strains,  and  stresses 
at  each  joint. 

Inelastic  Range 

The  procedures  described  in  the  elastic  range  may  be  extended  to  cover  the  inelastic 
range  by  introducing  the  variable  inelastic  modulus,  E',  when  the  stress  in  the  section 
exceeds  the  yield  stress.  Equation  B-16  becomes 


6.' 

) 


M  L 

=  "FT 

} 


(B-40) 


The  solutions  can  be  obtained  by  small  increments  of  stress  at  each  joint.  When  the  yield 
stress  is  reached,  the  inelastic  modulus  becomes  effective  at  that  joint.  The  forces, 
moment,  stresses,  strains,  and  deflections  are  calculated  at  the  point  by  the  equations 
shown  in  the  elastic  program.  The  stress  Increments  can  continue  until  the  ultimate  stress 
is  reached. 

For  each  change  of  slope  in  the  stress-strain  curve,  a  corresponding  yield  stress  is 
given.  Since  the  strain-hardening  materials  do  not  have  a  definite  plastic  range,  it  becomes 
necessary  to  represent  the  stress-strain  curve  by  a  series  of  lines  (more  than  two).  Figure 
B-8  shows  the  idealized  stress-strain  curve  for  the  circular  arch  tested.  The  variable 
inelastic  modulus,  E',  is  represented  by  Ei  and  E2  in  this  case. 


Figure  B-8.  Idealized  stress-strain  curve  for  static  analysis. 
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Computer  Programs 

Computer  programs  (B-l  and  B-2)  were  used  for  the  analysis  of  the  circular  arch  for 
the  elastic  and  Inelastic  ranges.  The  programs  follow  the  theoretical  derivation  of  equations 
given  in  this  appendix.  The  source  programs  are  written  in  FORTRAN  language  and  can  be 
carried  out  on  an  IBM  1620  or  IBM  1Q20II  computer.  A  flow  chart  showing  the  sequence  of 
operations  of  Program  B-l  is  given  in  Figure  B-9.  Program  B-2  is  relatively  simple,  and 
no  flow  chart  Is  needed.  The  Input  specifications  for  the  programs  are  as  follows: 


Number  of  bars: 

Stiffness  parameters: 
Elastic  support  constant: 
Loading: 


variable  but  subject  to  the  limitations  of  the  computer 
memory  system 

can  be  varied  from  joint  to  joint 
variable 

uniformly  distributed  over  any  portion  of  the  arch 


Presentation  of  Numerical  Results 


The  static  properties  for  materials  used  in  the  arches  are  given  in  Table  2.  Due  to  the 
strain-hardening  properties  of  some  specimens,  particularly  the  bottom  flanges,  the 
idealized  stress -strain  relationship  is  represented  by  the  trilinear  curve  of  Figure  B-8 
instead  of  by  the  usual  bilinear  stress-strain  curve. 

The  average  elastic  modulus  of  27  x  106  psi  was  used;  the  inelastic  moduli,  Ei  and 
£2,  were  0.4E  and  0. 08E,  respectively.  The  average  ultimate  stress  for  the  bottom  flange 
was  72,  500  psi.  Other  input  parameters  for  static  analysis  of  the  arch  were  as  follows: 


Number  of  bars 

n 

= 

16 

Angle  of  segments 

= 

0. 1  j82  radian 

Offset  angle 

0 

= 

0.06117  radian 

Radius  of  arch 

R 

= 

96. 0  in. 

Stiffness  parameter 

EI 

s 

407. 1  x  106  lb-in. 2 

Stiffness  parameter 

EA 

= 

177.  4  x  106  lb 

E  of  top  flange 

Eu 

- 

27. 9  x  106  psi 

E  of  bottom  flang: 

E! 

= 

26. 7  x  106  psi 

Initial  misalignment 

= 

0. 0195  in. 

Elastic  support  constant 

k 

- 

6. 25  x  10'6  in.  /lb 

Upper  fiber  from  neutral  axis 

cu 

= 

1.29  in. 

Lower  fiber  from  neutral  axis 

cl 

= 

2.71  in. 

Loading  was  considered  uniformly  distributed  over  one-half  the  arch  length.  The 
magnitudes  of  the  uniform  loads  used  In  the  analyses  were  increased  by  increments  of 
100  lb/in.  from  zero  to  500  lb/ in.  The  yield  stress  of  44, 000  psi  was  first  reached  in  the 
bottom  flange  at  the  1/4-point  of  the  arch,  where  the  stresses  caused  by  axial  force  and 
flexure  were  of  the  same  sign;  the  load  causing  this  yielding  was  507  lb/in.  At  a  load  of 
545  lb/ in. ,  the  yield  stress  was  developed  in  the  bottom  flange  at  the  3/4-point.  Higher 
loads  of  680  and  704  lb/in.  produced  stresses  of  58,000  psi  (point  B  in  Figure  B-8)  in 
the  bottom  flanges  at  the  1/4-point  and  3/4-point,  respectively. 
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START 


Figure  B-9.  Flow  chart  for  the  static  analysis  of  the  circular  arch. 
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r>  a  ri  r>  rt  r>  r\  r>  nnriA  n  a 


Program  B-l.  Analysis  of  the  Statically  Loaded  Arch,  Part  I 


ANALYSIS  OF  ARCHES  »Y  VIRTUAL  WORK  METHOD 
PROGRAMED  /  RICHARD  H,  CHIU  /  MARCH  1966 

DIMENSION  XIII?) *$Mf 17) «EWI 171 »X  I 17) ,  Y I 17) »XX ! 17) , YY ( 17) ,  W  t 17) 
DIMENSION  DRXI171 ,DRYi 1?)  ,FXM(  17)  * T  <1  ?  J  , ST  1 1 7 , 2  ( 

1  FORMAT  1IA«6F12,6«1A) 

2  FORMAT!! 3.5F15.A) 

3  FORMAT! I 3.5F15.8) 

*  FORMAT | /6H  J01NT7X1HA1AX1HX1 AX 1HY 1 AX2HXX 1 3X 2HYY/ 1 

5  FORMAT ( /6H  JO t NT7X2HXW1 3X2HY W 1 2 X 3HSXW 1 2 X3HS yw 1 2X2 HXM/ | 

6  FORMAT  (  /6H  JOINT6XAHE  ,rt.  10XAH  SFW12XAHDEFX/ ) 

7  FORMATI/6H  JOl  NT7  X 1 HH1  AX  lHV  1 2X  *>HEW  I  1  >  1 2X  3HSYM)  3X2HMX  /  J 

8  FORMAT  I6HJOINT  7X1HH t AX  1 HV I  AX 1HS 1  AX IhT 1 3XAHM ! 1) ) 

10  FORMAT I/6HJ0INT  6X2H0X1 3X2MD Y1 3X2HDR 1 3X?HDT 1 0X6HFXM I  I ) / ) 

102  FORMAT  !6HJ0INT  6X2HF T1 3X2HKB1 3X2HST 1 3X2HSb  ) 

100  READ  l,N*uA*&.R.CSA,£C.SPRC.Nl 

READ  1,  N3,DW»C1 ,C2,E»E1 
GEOMETRY  OF  CIRCULAR  ARCH 
00  9  1  *1  «N,6 

9  READ1 ,N2.X I <l>,XI(I  +  l>,XI(!+2!,xm  +  3),XI<!,A>,xm  +  S) 

00  110  1*1, N 

110  W! I 1*0. 

X!  1 1*0. 

Yi  1  »«0. 

XXII). 0. 

YY! 1)«0, 

EST ABL I SH  THE  COORDINATES  OF  JOINTS  1  THRU  N 

11  A*DA 
XM*0. 

PUNCH  A 

13  00  1A  1-2, N 

XII l-R*«COSIbl*!l.-CO$IA| )-SIN(A|*SIN(B)) 

Y(!  |-R*ISINI6)*SIN(A-Bn 

SECTIONAL  COMPONENTS  IN  X  AND  Y  DIRECTIONS 

15  XXI 1 )-X( l)-X( 1-1) 

YYIII-YI I)-Y(I-l) 

PUNCH  2,  1 ,A.X1 I ) ,YI I ) ,XX! I) ,YY( I | 

1 A  A-A+OA 

X  AND  Y  COMPONENTS  OF  THE  APPLIED  LCACS 
MOMENTS  ABOUT  0 
PUNCH  5 
DO  19  J* 1 ,N1 
SXW-O. 

SYW*0. 

DO  16  I-N3.N 
W( I )*W( I  I +0W 
YW-XXI I )*WI II 
XW* YY I  1 )*WI  I  I 
SXW-SXW-fXW 
SYW-SYW4YW 

XM«XM+XW*( Y! I  ) -YY 111/2*1 +YW* 1 X ( I 1 -XX! I  1/2.1 
C  PUNCH  2,  I  , XW , YW , SXW , SYW , XM 

16  CONTINUE 
PUNCH  2,  J , W I N ) 

VN-XM/X  <  N ( 

VO-SYW-VN 

HO* SXW 
RSX*0* 

RSY-O. 

CO  19  K* 1 . 3 
XM«0. 

SMX*0, 

SYM=0, 

A—  B 

00  20  1*2, N 
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Program  B-l.  Analysis  of  the  Statically  Loaded  Arch,  Part  I  (Contd) 


MOMENT  AT  SECTION  X 

ACCUMULATIVE  AREA  OF  THE  MOMENT  DIAGRAM 
OS*R*OA/(E*XI ( I  I ) 

1 F  <  K-2  )  118.117,117 

117  XM*XM*Y Y ( I l*(HO-W<  I  I *YY< I  1 /2. )*XX( I )«( VO-Wt  II*XXI|I/2,I 
GO  TO  119 

1 18  XM  «  Xt I >-Yt I > 

119  EW< I »*XM«DS 
SMX*SMX ♦EW f I ) 

SYM*SYM+EW( I)*X( I ) 

IFIK-2I  20,121.121 

121  HO«HO-W( I )*YY( I  I 
VO«VO-W(  I  ) *XX (  I  1 
IF  IK-2)  20.20,17 

17  A*A+DA 

ORXt I )«-HO*XX( I ) / ( CSA*E ) 

R5X  *  RSX  +  DRX ( I ) 

ORY ( I )«-V0*YY( I ) / ( CSA«E  » 

RSY«RSY+ORY( 1 ) 

RR«HO*COS ( A  I +VO*S I N ( A  I 
RT-VO*COS(A)-HO*SlN( A) 

T ( | ) *RT 

FXM ( I )*XM 

STC I.1)«RT/CSA  +  XM*C1/XI (  I  I 
ST ( I  ,2i*RT/C$A-XM*C2/XI  C  I  ) 

25  PUNCH  2.  I.HO.VO.RR.RT.XM 
GO  TO  20 

18  PUNCH  2,  I .HO.VO.EWI I ) .SYM.XM 
20  CONTINUE 

131  EWCN) *-SYM/X»N) 

EW( 1)«-I EW( Ni+SMX ) 

SXM»0. 

SEW«EW ( 1 1 
PUNCH  6 

PUNCH  3.  N2.EWI 1 ) , SEW.SXK 
00  24  I  * 2  »N 
SXM*SXM+SEW*YY( I ) 

C  PUNCH  3,  I  .Ea'I  I  )  .SEW.SXM 

2*  SEW*SEW+EWC I  I 

IFOC-21  22,21,141 
C  DEFLECTION  OF  THE  arch 

141  A=-R 

PUNCH  102 
DC  26  ; *2 ,N 
ST  T  «ST ( I ,1 >/E 
STU  =  0T ( I ,2  > /El 

26  PUNCH  3,  |  ,ST(  I  .1  1  ,r,T  (  I  ,2)  .STT.STb 
DX  =  0. 

Dy*0.  * 

SF  W*EW ( 11  f 

ROX  a-R  S  X / Y ( N I 

ROY  *— RS  Y / X  I N I 

PUNCH  10 

DO  2  3  I  *  2 , N 

DX*DX-,SFfc*YYIt)+ROX*TY{l!+ORXIl! 

DY  =  DY»SEW*XX ( t  f»ROY*XX«  I  l+ORYI  I ) 

SEWaSEW+Ew  I  I  ) 

Aza+da 

DR*DX*COS(  A),DY*S  |M  a  i 
DT  =DY*COS(A)~i'X*SIN(  A  I 
C  EFFECTS  OF  STRESS  AMPLIFICATION 

AMPF«FXM<  IJ/CEX-IIII  +  TII)  •  CDR-IC  I  ) 

FXM(  |  ) *  F  XM  <  I  J *AV.PF /  1  i.L  . 

IF  (AMPF-0.5)  142,142.143 
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Pr  ogram  B-l.  Analysis  of  the  Statically  Loaded  Arch,  Part  I  (Contd) 


14?  ADX  *  ox 
ADY*DY 
GO  TO  23 
143  AOX«DX*AMPF 
AOY=OY*AMPF 
DR*DR*AMPF 
OT»OT*AMPF 

23  PUNCH  3.  I  ,ADX.Af)Y.r>P*OT,FXMI  1) 

GO  TO  19 

EFFECTS  OF  HORIZONTAL  ELASTIC  SUPPORT 

21  XR*  ( SXM*3PRC/SXM1 — SXH1/SXM1 
HO=SXW-XR 

VO*  XR+VO+SYW 
PUNCH  8 

AXIAL  AND  SHF AR  STRFSSES  IN  THE  ARCH 
133  RR=HO*COS <  B ) -VO*3 I N  (  8  } 

R T  *HO* $  I N ( b ) ♦ VO*C OS ( S  > 

XM*0  • 

PUNCH  <.  *  N2  «H0  « VO  *RR  *  RT  »XM 
GO  TO  19 

22  SXM1=SXM 
PUNCH  7 

19  CONTINUF 

1 F ( SENSF  SWITCH  1)  100,101 
]01  PAUSE 
END 


Program  B-2.  Analysis  of  the  Statically  Loaded  Arch,  Part  n 


DIMENSION  X(9,17).SX<9.9),R(4),5R(4I 

1  FORMAT  C3X.4F15.6) 

2  FORMAT  (33X.3F15.6l 

3  FORMAT  (  I4,2F8.0.3F9.0.2F8.0,2F8.4) 

4  FORMAT  (3X.4F15.8/I 

5  FORMAT  <  //3X| 

6  FORMAT  ( /AH  BAR,5X.1HW.6X,4HR( 11 .4X.4HR12) »5X.4HR(3) .5X,4HR(4) , 
14X.5HV($R)/l 

7  FORMAT  (/4H  JT..2X.6HDIR.  F . 3X . 5HSHE AR « 3X . 6HMOMENT . 7X . 8HSTRESSES « 
llOX,7HSTRAtNS.t>X,  lortDEFLECT  ION) 

6  FORMAT  {  4M  NO. ,3X,4M(L6 I ,5X.4H(Lbl ,9H  (  IN-Lfll  .3X.6HT  <  PS  1  1  ,3X, 

16HB(PSt).17H  T  IN/IN  b  IN/ IN.2X.5HRC IN) ,3x»5Hl t IN)/ > 

9  FORMAT  ( 3X.  F15.8/ I 
DO  111  J*1.9 
DO  111  1-1.9 

111  SX< I . J|«G. 

DO  112  1*1.4 

112  SR ( I  1=0. 

sw*o. 

100  READ  3,  M.C5.C13.S 
READ  9  , W 

READ  1 ,91 3) »R(4) ,X(2» 11 .XI 1 .11 
DO  11  J  =  2 ,16 

1 1  READ  2,X(2.JI.X(1«J| 

READ  4,R( 11 ,R< 2) .X( 2 . 17) ,X< 1 ,17  I 
V*X<  2,1 l+0.763*( X( 2 ,2 )-X<2 .1  I  1 
DO  12  1*1,4 

12  R( 1  I =A8S  <  R ( I  I  I 
DO  21  J=2.4 

21  READ  1,  X(4.J),X(  5.J) ,X(6.J) ,X(7,J) 

READ  1,  X(4,5i ,X( 5.5) ,X(6«5l .XI7.5I 

X(4,5)*X(4,5)*C5 

X(5,5)=X(5,5)*C5 

DO  22  J  =  6 » 1 2 

22  READ  1,  X(4,J)  , X ( 5  ♦  J  )  ,X(6«J)  «  X  (  7  «  J  ) 

READ  1,  X(4,13),X(5,13),X(6.13),X(7,13) 

X(5,13)=X(5.13)*C13 

X(4,13)*X(4,13)*C13 

DO  23  J=14,17 

23  READ  1.  XI4.J) ,X(  5.JI .XI6.JI ,X(7,J| 

DO  24  1=4,7 

DO  25  J  =  2 ,17 
25  X( I , J I =  -X ( I ,Ji 
2A  XII ,1 !=X( 1,17) 

READ  8 
DO  31  1*8,9 

31  X ( I ,1 ) *0. 

X ( 3 , 1 ) *0 , 

DO  32  J  =  2  » 1 7 

32  READ  2,  XI8.J) »X I  9, J ) , X I  3  ,  J ) 

DO  42  J* 1 , 1 7 

DO  41  1=6.7 

41  X(l.J)*X(!«Jt*i.C£+6 

42  X(3,J)=XC3,J1*!.0E*3 
GO  TO  (45, 31. 52). M 

45  Sw  =  w 
SV*V 

DO  44  1*1,4 
44  SRC !)*R(  I  ) 

DO  43  J=1 ,9 
DO  43  1*1,9 
K*l+2«( J-l ) 

43  SX (  I  ,J)=X<  I  ,X) 

GO  TO  60 
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Program  B-2.  Aialysis  of  the  Statically  Loaded  Arch,  Part  n  (Contd) 


51  FAC«-(S-*SX(5.3l  >/X(5»5> 

GO  TO  53 

52  FAC-  (  S-SX  (5,7  J1/XC5.13) 

53  W-FAC-W 
V-FAC-V 

00  5*  1-1.4 

54  Rm>R(I  ) *FAC 
DO  55  J-1.17 
00  55  1-1,9 

55  XU*J)«X(I,J)*FAC 
SW-SW+W 
SV-SV+V 

DO  57  1-1,4 
57  SR< Il-SRI I ) +R ( I ) 

00  56  J-1,9 
00  56  1-1.9 
X«H-2*CJ-ll 

56  SX( I « J) -SX ( 1 ,  J 1 +X ( I ,X  j 

60  NJ- <  J-2 ) *2 
PUNCH  6 

PUNCH  3,  NJ,  W,  R ( 1 ) ,  R ( 2 ) ,  R(3),  R(4J,V 

PUNCH  7 

PUNCH  8 

00  61  J-1,17 

X-J-l 

61  PUNCH  3,  K,  X (  1  •  J )  •  X  <  2  ,  J  )  ,  XI3.J),  X(4,J),  X ( 5  ,  J  )  ,  X(6,J), 

1  X ( 7, J) .  X ( 8, J I ,  X ( 9 , J ) 

PUNCH  6 
NJ-NJ/2 

PUNCH  3.  NJ.SW.SR < 1 ) ,SR( 2l ,SR( 3) »SRI4)  ,SV 

PUNCH  7 

PUNCH  8 

DO  62  J-1,9 

X-J-l 

62  PUNCH  3,  X  »SX ( 1 , J ) »SX  (  2  »  J  ) »SX  (  3.J) »SX  ( 4.J) . £X ( 5 . J I ,SX ( 6 , J I . 
1SX(7.J) ,SXI B.J) ,SX(9«J) 

IF  (SENSE  SWITCH  2)  100,101 
101  CALL  EXIT 
END 
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A  load  of  866  lb/in.  developed  the  ultimate  stress  of  72,  500  psi  in  the  bottom  flange 
at  the  1/4-point.  At  a  load  of  877  lb  in. ,  the  bottom  flange  at  the  3  4 -point  also  attained 
the  ultimate  stress.  If  this  load  were  maintained,  the  arch  would  continue  to  deflect  until 
collapse  occurred. 

The  results  of  the  theoretical  response  in  the  elastic  and  inelastic  ranges  are  given 
in  Tables  B-l  through  B-4.  Table  B-l  summarizes  the  reactions  and  radial  shear  at  strain 
rosetto  3^-13-14-15.  The  location  of  the  rosette  is  shown  in  Figure  6;  the  directions  of  the 
reactions  are  shown  in  Figure  B-l. 

Table  B-2  summarizes  the  theoretical  results  at  the  1/4-point  of  the  arch.  At  this 
section  the  direct  forces  are  maximum  and  the  shearing  forces  are  minimum  for  each  level 
of  load.  Maximum  stress  and  strain  values  occur  at  this  location  in  the  bottom  flange. 
Positive  radial  deflections  are  in  the  outward  direction. 

Table  B-3  summarizes  the  theoretical  results  at  the  c-own  of  the  arch.  The  positive 
shear  and  tangential  displacements  are  maximum  at  this  location.  The  direction  of  positive 
tangential  displacement  is  shown  in  Figure  B-3. 

Table  B-4  summarizes  the  theoretical  results  at  the  3/4-point  of  the  arch.  The  bending 
moments  and  the  radial  deflections  are  maximum  at  this  location;  the  axial  and  shearing 
forces  are  small. 


Table  B-l.  Summary  of  Elastic  and  Inelastic  Reactions  and 
Shears  of  the  Statically  Loaded  Arch 


Load 
(lb,/ in. ) 

R(l) 

(lb) 

R(2) 

(lb) 

R(3) 

(lb) 

R(4) 

(lb) 

V(SR)v 

(lb) 

Elastic 

100 

4,590 

2,520 

1,800 

4,860 

-1,700 

200 

9, 180 

5,040 

3,600 

9,720 

-3,400 

300 

13,800 

7,560 

5,400 

13,580 

-5,100 

400 

18,400 

10, 080 

7, 200 

19,440 

-6,800 

500 

23,000 

12,600 

9, 000 

24,300 

-8,500 

507 

23, 300 

12,800 

_ 

9,120 

_ 

24,600 

-8,630 

Inelastic 

545 

25, 000 

13,800 

9,760 

26,400 

-9,  230 

680 

31, 200 

17,200 

12,  200 

33,000 

-11,600 

704 

32, 200 

17,900 

12,500 

34,100 

-11,900 

866 

39, 700 

21,900 

15,500 

42,  000 

-14,700 

877 

40, 200 

22,200 

15,700 

42,600 

-14,  900 

*V(SR)  are  shearing  forces  at  strain  rosette  SG-13- 14-15. 
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Summary  of  Elastic  and  Inelastic  Responses  at  the  Crown  of  the  Statically  Loaded  Arch 
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Computations  of  direct  force,  shear,  bending  moment,  stresses,  strains,  and 
deflections  were  made  for  all  joints.  However,  Tables  B-2  through  B-4  tabulate  those 
values  for  only  the  1  4-point,  the  crown,  and  the  3  4- point.  The  distribution  of  direct  forces 
along  the  arch  varied  from  a  maximum  at  the  1  4 -point  to  a  minimum  at  the  3  4-point  (e.  g. , 
from  10,  400  lb  to  8,  500  lb  at  a  load  of  200  lb  in. ).  The  1  4-point  and  3  4-point  had  maxi¬ 
mum  bending  moments  and  nearly  zero  shear.  The  fact  that  the  shear  was  not  zero  at  these 
points  indicates  that  the  maximum  bending  moment  does  not  occur  exactly  at  either  quarter 
point.  This  would  tend  to  provide  a  theoretical  stiffness  of  the  arch  that  is  slightly  higher 
than  it  would  be  if  a  greater  number  of  bars  was  considered  in  the  theoretical  solution. 

A  theoretical  curve  of  load  versus  deflection  is  shown  in  Figure  23. 
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Appendix  C 

DYNAMIC  ANALYSIS  OF  A  CIRCULAR  ARCH 


The  dynamic  response  of  a  two- hinged  circular  arch  subjected  to  dynamic  load 
distributed  uniformly  over  one-half  the  arch  length  was  calculated  for  the  elastic  and 
inelastic  ranges.  A  substitute  discrete  framework  was  used  in  the  analysis,  and  the 
numerical  results  were  obtained  by  the  iterative  method  of  numerical  integration.  The 
original  investigation  was  performed  under  contract  by  Eppink  (References  15,  16,  and  1?). 
Variations  and  additional  comments  were  made  primarily  for  clarity  and  continuity  of  this 
report. 

Analogous  Framework 

As  pointed  out  in  the  static  analysis  of  the  circular  arch  of  Figure  B-l,  the  exact 
solution  of  the  continuous  arch  with  distributed  mass  is  impracticable.  This  is  even  more 
so  in  the  dynamic  analysis  of  the  circular  arch.  In  order  to  simplify  the  solution,  it  is 
desirable  to  represent  the  original  continuous  arch  by  a  discrete  framework  consisting  of 
a  series  of  rigid  bars  and  flexible  joints.  The  following  assumptions  are  made  regarding 
the  analogous  framework  in  the  elastic  range: 

1.  The  rigid  bars  are  assumed  to  be  massless  and  straight. 

2.  The  flexible  joints  have  the  elastic  properties  of  the  arch  section  included  between 
the  midpoints  of  the  adjoining  segments. 

3.  The  distributed  mass  of  the  segment  is  lumped  into  a  concentrated  mass  at  the 

joint. 

Due  to  the  complexity  of  the  stress-strain  relationship  in  the  inelastic  range,  further 
assumptions  regarding  the  section  properties  of  the  arch  are  needed: 

1.  The  cross-sectional  area  of  the  arch  is  composed  of  two  flanges  and  a  thin  web  of 
zero  area. 

2.  The  web  is  assumed  to  be  rigid  in  shear  but  resists  no  axial  force. 

Schematic  representations  of  the  circular  arch  and  the  analogous  framework  are  shown 
in  Figures  C-l  and  C-2,  respectively.  The  joints  of  the  discrete  framework  are  numbered 
consecutively  from  0  at  the  left  support  to  n  at  the  right  support,  and  the  bars  are  numbered 
from  1  to  n.  A  subsc  ript  notation  is  used  to  identify  the  locations  of  joints  along  the  arch. 

A  representative  joint  is  designated  by  subscript  j;  joints  j-2,  j-1,  j+1,  j+2,  and  j+3  are 
designated  by  h,  i,  k,  I,  and  m,  respectively.  The  same  system  is  also  used  for  bars  such 
that  bar  j  connects  joints  i  and  j. 

The  deformed  arch  is  represented  by  the  radial  and  tangential  components,  djr  and 
dj(-,  of  the  displacement  of  the  joints.  All  displacements  are  measured  from  the  undeformed 
position  of  the  joints,  and  the  directions  of  positive  displacements  are  shown  in  Figure  C-2. 
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Development  of  Equations  of  Motion 

The  forces  acting  at  the  ends  of  a  typical  bar  and  Joint  are  shown  in  Figure  C  >3.  They 
Include  the  axial  force,  ,Pjt;  the  transverse  shearing  force,  P<r;  the  bending  moment,  Mj; 
the  damping  forces,  cedjt  and  ceijr;  and  the  applied  force,  Pj.  The  applied  force  is 
expressed  in  terms  of  its  components,  P,1  and  Pjr,  to  the  left  and  right  of  the  joint  j, 
respectively.  Summing  the  components  of  all  internal  forces  in  the  radial  direction,  one 
obtains  the  following  expression  for  internal  resulting  force  in  the  radial  direction,  Zjr 

Zjr  =  -  Pjt  sln  (I  '  *i)  -  Pkt  sln  (I  +  *k) 

+  Pjr  cos  (i  -  *j)  •  Pkr  cos  (l+  *k) 


where  denotes  the  rotation  of  bar  j  and  ^  denotes  the  rotation  of  bar  k.  Considering 
a  small  deformation  such  that  sin  \J>  »  ^  and  cos  ^  a  1,  the  above  equation  can  be  written 
in  the  following  form: 


V  ■  •(%  *  pw)s‘"f  *  (Vj  - 

*  (PJr  •  Pkr)  C0B  I  *  tPjr  *  Pkr*k)slnf 


(C-l) 


Similarly,  the  tangential  component  of  the  internal  resulting  force,  Zjt,  may  be 
obuined  by  summing  all  the  internal  forces  in  the  tangential  direction  as  follows: 


Zjt  =  -Pjtcos 


(I  '  *j)  +  Pktcos(l  +  *k) 


-  Pjrsln\2 


(I'*  j)'  Pkr8in(i  +  *k) 


and  for  a  small  deformation  this  reduces  to 


z„  ■  -(p)t  -  pkt)"»l  -  (pit*i  ♦  pkt'sk)s"'l 

—  (p  *  P  \  Sin  ^  /p  it  -  P  ■/,  \  r»ne 

\*jr  *  kr/  sin  2  V *jr  vj  *  kr  v  k/  2 
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Z2r 


Figure  C-3.  Free  body  diagrams  for  a  typical  bar  and  joint. 


The  shearing  forces  Pjr  and  Pgr  in  the  above  equations  are  related  to  the  bending 
moments  Mi,  Mj,  and  Mfc  as  follows: 


where  6\  and  are  deformations  of  the  bar  j  and  bar  k,  respectively.  The  axial  forces, 

St,  ana  shearing  forces,  Pjr,  are  taken  at  the  center  of  the  bars,  and  the  bending  moments, 
j,  are  taken  at  the  centers  of  the  joints. 
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The  external  pressure  is  considered  to  act  normal  to  the  arch  surface.  The  distributed 
load  is  represented  by  a  series  of  concentrated  forces  acting  at  the  joints.  For  convenience, 
the  concentrated  forces  are  represented  by  their  components,  Pj*  and  Pjr,  which  are 
positive  in  the  outward  direction.  By  summing  the  external  forces  in  the  radial  direction, 
one  obtains  the  expression  for  the  external  resultant  force  in  the  radial  direction,  Wjr, 
as  follows: 


Wjr  =  P/  cos  (I  '  *])  +  Pjr  cos  (4  +  *k)  -  ce  V 


where  ce  is  the  coefficient  of  external  viscous  damping,  and  a  dot  superscript  denotes 
differentiation  with  respect  to  time.  In  the  absence  of  specific  information  regarding  the 
internal  damping  of  various  types,  the  external  damping  is  taken  to  be  the  total  damping  of 
the  system.  The  coefficient  ce  is  obtained  from  experimental  sources.  For  small  displace¬ 
ment,  the  above  equation  becomes 


Wjr  '  (P/  *  Pir)  “8  $  *  (Pi‘  *)  -  Pi'  *k) Bl"  i  -  V  <c"» 


Similarly,  by  summing  the  external  forces  in  the  tangential  direction,  one  obtains 
the  expression  for  the  external  resultant  force  in  the  tangential  direction,  Wjt,  as  follows: 

W.t  =  -P*  sin  (|  -  *.)  ♦  P/sind  +  *k)  -  cedjt 


and  for  a  small  displacement,  it  becomes 


where  iijr  and  dijt  are  the  radial  and  tangential  accelerations  of  the  mass  mj. 

The  schematic  representation  of  the  displacement-deformation  relationships  is  shown 
in  Figure  C-4,  where  ab  and  a'b'  are  the  positions  of  bar  j  before  and  after  deformation. 

Let  the  projection  a’c'  on  ab  be  designated  by  L  +  dy  then  dj  can  be  expressed  in  terms 
of  the  displacements  of  joints  i  and  j  as  follows: 

Jj  =  (4jt  "  4it)cosl  +  (4jr  +  *ir)ain  4  (C-8) 
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Figure  C-4.  Schematic  representation  of  displace  ment-deformation  relationships. 

Let  the  projection  c'b'  in  the  direction  normal  to  ab  be  designated  by  L^j;  then 

*j  =  r[(4jt  +  Ait)sini  -  (4jr  -  4ir)c08f]  (C'9) 

The  length  of  the  deformed  bar,  a’b',  can  be  expressed  in  terms  of  the  deformation,  ij, 
as  follows: 


*  ‘  L  V*i2  *  (‘  ’4 


The  deformation  of  the  bar  j  is  then 


=  L 


-1 


* 


(C-10) 


Similarly,  the  rotation  of  the  bar  j,  ^ j,  can  be  expressed  in  terms  of  Jj  and  0 j  as  follows: 


\}> .  =  tan 


1  +-T 


(C-ll) 
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Similarly, 


sin  ^  = 


*4 


V*i 2*  ('*4)' 


*1  (>  -  4) 


cos  ^  = 


*4 

1  +  L 


(‘  *4) 


* 1  -  iV 


where  ij  and  are  the  first-order  approximations  of  ij  and  0j.  The  angular  change  at 
joint  j,  £j,  is 


9.  -  \f>.  - 
j  *k 


(C-12) 


and  its  first-order  approximation,  flj,  is 


9.  =  -  4>. 

]  **  J 


(C-13) 


The  geometric  relationships  and  the  equations  of  motion  developed  in  the  preceding 
paragraphs  are  general;  they  are  valid  irrespective  of  the  material  properties.  In  order  to 
formulate  the  force-deformation  relationship,  distinctions  between  the  material  properties 
in  the  elastic  and  inelastic  ranges  must  be  made. 

Elastic  Range.  Taken  from  midpoints  of  the  adjacent  bars,  the  deformation  to  the 
left  and  to  the  right  of  joint  j  can  be  expressed  in  terms  of  internal  forces  and  section 
properties  of  the  joint  as  follows: 


6. 


1  LPit 


j  2  E  A 


6r 

i 


1  LPkt 


2  EA 


‘j 


(C- 14) 


where  E  Aj,  the  extensional  stiffness  parameter,  is  the  product  of  the  elastic  modulus  and 
cross-sectional  area  of  joint  j.  The  total  deformation  of  bar  j  is  expressed  as  follows: 


'i  ■  V  *  V 


(C-15) 
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Combining  the  above  expression  with  Equation  C-14  and  solving  for  the  axial  force  due 
to  elastic  strain,  Nj,  one  obtains 


In  the  absence  of  extensional  damping,  the  total  axial  force,  Pjt,  equals  Nj. 

Similarly,  the  bending  moment,  Mj,  is  expressed  in  terms  of  angular  change,  Sj, 
and  the  flexural  stiffness  parameter,  Elj,  of  joint  j  as  follows: 


(C-17) 


Inelastic  Range.  In  the  inelastic  range,  the  actual  arch  cross  sections  are 
approximated  by  a  section  consisting  of  two  flanges.  For  simplicity,  the  idealized  bilinear 
stress-strain  curve  of  Figure  C-5  is  used  to  correlate  the  force-deformation  relationships. 
The  curve  is  shown  in  dimensionless  form  such  that  the  slope  of  the  elastic  portion  is  unity 
and  yielding  occurs  when  the  value  of  f'fy  is  equal  to  one.  The  slope  of  the  inelastic  portion 
is  E}  E,  the  ratio  of  the  inelastic  modulus  to  the  elastic  modulus  of  the  arch  material. 
Unloading  is  assumed  to  take  place  along  a  line,  A  B,  parallel  to  the  elastic  line.  The 
intercept  of  A  B  on  the  e  ey  axis  is  the  instantaneous  set,  b/ey.  Unloading  will  continue 
along  this  line  until  a  new  "elastic"  limit  is  reached  at  B.  The  projected  length  of  A  B  on 
the  vertical  axis  is  2f/fy.  Further  unloading  will  occur  along  a  line  with  a  slope  of  E\,E. 

A  possible  path  of  a  loading  and  unloading  cycle  is  represented  by  OABCD  in  Figure  C-5. 


Figure  C-5.  Bilinear  stress-deformation  relationship. 
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In  the  inelastic  range,  the  stress-deformation  relationship  is  rather  complicated.  It 
becomes  necessary  to  determine  the  deformations  in  the  flanges  before  the  stress,  moment, 
and  axial  force  in  each  joint  can  be  calculated.  A  schematic  representation  of  the  flange 
deformations  is  shown  in  Figure  C-6.  Additional  subscript  notations,  t  and  b,  are  intro¬ 
duced  to  designate  the  top  and  bottom  flanges,  respectively. 


Figure  C-6.  Schematic  representation  of  deformed  joint. 


To  avoid  the  time-consuming  process  of  solving  numerous  simultaneous  equations, 
additional  approximations  are  needed.  The  distribution  of  angular  change  to  the  left  and  to 
the  right  of  a  joint  is  inversely  proportional  to  the  relative  stiffness,  EI/L,  of  the  adjacent 
bars;  the  distribution  of  axial  deformation  is  inversely  proportional  to  the  equivalent  cross- 
sectional  areas,  EA,  of  two  joints.  Considering  the  arch  to  have  uniform  cross  section 
along  its  length,  the  deformations  are  given  by  the  following  equations: 


•i‘  ’  l(‘j  *  Vi) 

ajb*  '  t('i  '  cjb  ’’)) 

V  ■  i(ak  *  cit  *j) 

aibr  -  j(‘k  -  b)b «)) 


(C-18) 
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where  cjt  and  c<b  are  the  distances  from  the  neutral  axis  of  joint  j  to  the  top  and  bottom 
flanges,  respectively.  Equations  C-18  give  exact  results  for  all  deformations  in  the  elastic 
range.  Specialized  conditions  must  be  used  to  account  for  the  boundary  condition  at  the 
support.  For  hinged  supports  the  elastic  angular  change,  B,  is  zero.  The  following 
expressions  are  obtained  for  joints  o  and  n: 


(C-19) 


As  a  result  of  the  assumptions  used  to  obtain  Equations  C-18,  the  equilibrium 
conditions  at  the  joints  can  only  be  met  by  introducing  the  average  axial  forces  and  moments 
as  follows: 


pJt  - 1  (V  *  pibr  *  V  *  V) 


(C -20) 


=  2  fob*  Cjb  ‘  Pjt*  Cjt  +  Pjb"  Cjb  '  Cjt) 


(C-21) 


Method  of  Numerical  Integration 

Equations  C-6  and  C-7  may  be  expressed  in  terms  of  the  radial  and  tangential 
components  of  displacement  for  joints  1  through  n-1  by  using  the  relationships  established 
in  the  preceding  section.  The  solution  of  the  dynamic  response  of  the  arch  is  obtained  by 
satisfying  a  set  of  2  (n-1)  simultaneous  nonlinear  differential  equations  of  motion.  To 
accomplish  this,  a  method  of  iteration  and  numerical  integration  is  used. 

The  process  of  iteration  begins  at  time,  tn,  with  known  values  of  displacements, 
velocities,  and  accelerations  of  all  joints  in  the  system  and  proceeds  to  solve  for  the  corre¬ 
sponding  quantities  for  a  time,  tn+i,  at  a  small  time  interval,  dt,  from  tn.  It  should  be 
noted  that  the  initial  values  of  displacements,  velocities,  and  accelerations  may  be  estimated 
or  assumed.  However,  good  estimates  will  ensure  a  speedier  convergence.  A  step-by-step 
operation  of  this  method  will  be  shown  in  the  following  paragraphs.  Each  of  the  following 
operations  is  to  be  repeated  for  the  components  of  displacement  at  each  joint  of  the  system. 

1.  Assume  or  estimate  a  set  of  values  of  displacements,  velocities,  and  accelerations 
at  a  chosen  time  tn.  For  example,  at  initial  time,  t0,  all  quantities  may  be  taken  to  be 
zero;  at  time  tn  the  quantities  may  be  assumed  as  the  corresponding  quantities  at  time  tn-1. 

2.  Assume  a  function  for  accelerations  for  the  small  time  interval  and  evaluate  the 
values  cf  velocities  and  displacements  at  the  end  of  the  time  interval  in  terms  of  the  accelera¬ 
tions,  velocities,  and  displacements  at  the  beginning  of  the  interval.  To  accomplish  this, 

the  following  equations  obtained  from  Newmark  **  are  used: 
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|  where  x  may  be  taken  to  be  either  component  4jr  or  djt  of  displacement  at  joint  j.  The 

dots  denote  differentiation  with  respect  to  time;  the  subscripts  n  and  n  + 1  correspond  to 
times  tn  and  tn  + 1,  respectively;  0  depends  on  the  assumed  variation  of  the  accelerations 
within  the  time  interval,  dt.  0  Is  1/6  for  a  linear  -ariation  of  the  accelerations. 

3.  Calculate  4j,  if/j,  and  of  Equations  C-10,  C-ll,  and  C-13  using  the  values  of 
displacement  at  the  end  of  the  time  interval. 

4.  Evaluate  Wjr  and  Wjt  of  Equations  C-4  and  C-5  for  external  force  components. 

5.  Evaluate  P«  and  M<  of  Equations  C-16  and  C-17  for  the  elastic  response  of  axial 
force  and  moment  at  joint  j.  The  corresponding  quantities  can  be  evaluated  from  Equations 
C-20  and  C-21  using  information  from  the  stress-strain  curve  and  Equations  C-18  and  C-19. 

6.  Evaluate  Pir,  Zjr,  and  Zj t  of  Equations  C-3,  C-l,  and  C-2  for  shearing  force 
and  the  components  ox  internal  forces  at  joint  j. 

7.  Determine  a  new  set  of  accelerations  of  Equations  C-6  and  C-7  with  the  results  of 
wJr>  wjt»  zjr.  and  zjt- 

8.  Compare  the  derived  accelerations  with  the  assumed  accelerations. 

9.  Substitute  the  derived  accelerations  for  the  corresponding  assumed  accelerations 
if  the  agreement  is  not  satisfactory,  and  repeat  steps  2  through  9.  Proceed  to  the  next  time 
interval  when  a  satisfactory  degree  of  agreement  is  reached,  and  repeat  steps  1  through  9. 

The  above  iterative  procedure  was  intended  for  the  dynamic  arch  response.  However, 
a  quasi-static  response  can  also  be  obtained  from  the  dynamic  program  by  introducing  a 
relatively  high  level  of  damping. 

Computer  Program 

A  computer  program  was  prepared  for  the  procedures  described  in  the  previous 
section.  The  source  program  is  written  in  ALGOL  language  and  can  be  carried  out  on  the 
Burroughs  B-5000  computer  at  the  University  of  Virginia.  The  general  i  labilities  and 
features  of  the  programs  are  as  follows: 

Number  of  bars:  variable  but  subject  to  the  limitations  of  the  computer 

memory  system 

Stiffness  parameters:  can  be  varied  from  joint  to  joint 

Support  conditions:  varying  from  fully  hinged  to  fully  fixed  supports* 

Loading  conditions:  arbitrary  in  distribution  and  magnitude.  The  time 

variation  of  load  is  considered  linear 

Output  quantities:  axial  force,  bending  moment,  radial  and  tangential 

components  of  displacement,  stress,  and  strain 

Properties  of  the  Circular  Arch  and  Load  Function 

The  geometry  of  the  two-hinged  circular  arch  and  the  distribution  of  loads  are  shown 
in  Figure  B-l.  The  parameters  of  the  arch  are  as  follows: 

Radius  of  arch  R  =  96. 0  inches 

Span  of  arch  L  =  143. 8  inches 
o 

Angle  of  opening  <t>Q  =  97. 01  degrees 
*See  Reference  15  for  derivations  of  moment- resisting  supports. 
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The  curvature  of  the  circular  arches  was  obtained  by  the  cold-roll  process,  and  the 
flanges  of  the  4M13  rolled  section  were  strained  beyond  the  yield  point.  The  dynamic  stress- 
strain  curves  do  not  show  the  familiar  yield  plateau.  For  simplicity,  the  Idealized  dynamic 
stress-strain  curves  for  the  component  parts  of  the  arch  cross  section  are  taken  to  be 
bilinear  as  shown  in  Figure  C-7. 

The  idealized  time  variation  of  the  dynamic  load  is  represented  by  a  triangular  decaying 
function  with  initial  peak  load,  p0,  as  shown  in  Figure  C-8.  The  load  duration,  t0,  is 
taken  to  be  1.6  seconds  for  all  magnitudes  of  the  peak  loads. 


Figure  C-7.  Idealized  dynamic  stress-strain  relationship. 
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Figure  C-8.  Load- time  relationship. 


The  cross-sectional  properties  of  the  arch  in  the  elastic  range  are  represented  by  the 
stiffness  parameters,  EA  and  El,  as  follows: 

EA  =  177.4  x  106  lb 

El  =  407.1  x  10®  lb- in. 2 


The  cross-sectional  properties  of  the  arch  in  the  inelastic  range  are  more  complicated; 
the  actual  cross  section  is  represented  by  two  flanges.  Attempts  were  made  to  provide  the 
closest  approximations  of  the  actual  cross-sectional  properties.  The  following  set  of 
parameters  were  used: 

3.  815  in. 2 
2.  288  in. 2 
30  x  106  pst 
28  x  106  psi 
55,  300  psi 
86,  500  psi 
3. 15  in. 

In  the  absence  of  specific  information  on  the  damping  characteristics  of  the  arch 
material,  it  was  assumed  that  the  external  damping  measured  by  experimental  means  repre¬ 
sented  the  total  damping  of  the  system.  An  average  damping  factor,  v,  of  0. 095  was 
observed  from  the  fundamental  antlsymmetrical  mode  of  the  free  vibration  tests.  In  order 
to  conform  with  the  damping  term  used  in  the  computer  program,  the  above  damping  factor 
was  converted  to  that  of  a  complete  ring  oscillating  at  the  frequency  corresponding  to  the 
"breathing"  mode. 


Area  of  top  flange 

Area  of  bottom  flange 

Modulus  of  elasticity  of  top  flange 

Modulus  of  elasticity  of  bottom  flange 

Yield  strength  of  top  flange 

Yield  strength  of  bottom  flange 

Distance  between  flanges 


At  = 
\  " 
Et  = 
Eb  = 

fyt  = 

£yb 
d  = 


/ 
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As  shown  in  Appendix  D,  the  first  antisymmetrlcal  frequency  of  the  arch  was  59.  9  cps, 
and  the  "breathing"  mode  frequency  of  a  complete  ring  was  317.  3  cps.  The  equivalent 
damping  factor,  v',  was  as  follows: 


v •  =  (0.095)  =  0.018 


Presentation  of  Numerical  Resu)*' 

Dynamic  responses  of  the  arch  in  the  elastic  range  were  obtained  for  peak  loads  of 
200  lb/in.  and  400  lb'in.  For  loads  above  400  lb/ in. ,  the  inelastic  behavior  of  the  arch  was 
anticipated  and  the  inelastic  program  was  used.  As  a  comparison,  the  load  of  400  lb  in. 
was  repeated  in  the  inelastic  solution.  The  loading  was  increased  in  increments  of  100  lb  in. 
up  to  700  lb/in.  and  then  increased  in  50-lb  'in.  increments.  A  maximum  dynamic  load  of 
800  lb/in.  was  used  in  the  solutions. 

The  results  of  dynamic  solutions  for  elastic  response  for  peak  loads  of  200  lb/in.  and 
400  lb/in.  are  given  in  Table  C-l,  along  with  the  maximum  values  of  bending  moments, 
axial  forces,  radial  and  tangential  displacements,  and  the  corresponding  times  when  the 
maximum  values  occurred.  The  inelastic  responses  for  loading  between  400  lb.  in.  and 
800  lb/in.  are  given  in  Tables  C-2,  C-3,  and  C-4.  All  solutions  were  based  on  a  16-bar 
system;  however,  data  are  presented  for  intervals  of  one-eighth  of  the  arch  length.  At  the 
maximum  load  of  800  lb  in. ,  the  convergence  of  radial  and  tangential  displacements  was 
slow.  The  maximum  displacements  were  not  obtained  at  several  locations  of  the  arch  as 
indicated  in  Table  C-4. 

The  maximum  elastic  and  inelastic  values  of  moments,  axial  forces,  and  radial  and 
tangential  displacements  at  the  quarter  points  and  the  crown  of  the  arch  are  presented  in 
graphic  form  in  Figures  C-9  through  C-12.  The  bending  moment  curves  for  the  quarter 
points  coincide  with  each  other  in  the  elastic  range,  but  their  signs  are  different.  Bending 
moments  at  the  crown  are  comparatively  low.  The  maximum  values  of  axial  forces  do  not 
vary  appreciably  at  the  quarter  points  and  at  the  crown  of  the  arch,  as  shown  in  Figure  C-10. 
The  maximum  radial  displacements  occur  at  the  3/4-point  and  the  maximum  tangential 
displacements  occur  at  the  crown  of  the  arch. 

The  time  interval  within  which  the  maximum  values  of  various  response  functions  are 
attained  is  only  a  very  small  fraction  of  th<  entire  load  duration.  The  decrease  of  the  load 
from  its  peak  value  is  very  small  within  this  time  interval. 

The  time-variant  traces  of  moment  and  displacements  show  that  the  major  contribution 
came  from  the  first  antisymmetrical  mode  of  vibration.  This  observation  can  be  made  from 
the  moment  curves  and  displacement  curves  of  Figures  C-13  and  C-14  for  a  peak  load  of 
200  lb  in.  The  period  of  major  oscillations  corresponds  to  the  period  of  the  first  antisymmet¬ 
rical  mode  of  the  arch  in  free  vibration.  It  should  be  noted  that  the  influence  of  the  third 
symmetrical  mode  is  also  present,  but  the  effects  are  comparatively  small. 

The  predominant  oscillation  of  the  axial  forces  oc-  <rs  at  a  frequency  corresponding  to 
the  second  symmetrical  mode  of  vibration,  which  is  close  to  the  "breathing"  mode  for  a 
complete  ring.  The  time  variant  traces  of  axial  forces  near  the  crown  and  supports  of  the 
arch  are  shown  in  Figure  C-15.  It  can  be  seen  that  the  variation  of  axial  forces  along  the 
arch  is  small. 

In  the  inelastic  range,  the  response  frequencies  and  modes  of  various  functions  are 
modified.  Increases  in  the  natural  period  of  the  response  curves  are  expected  as  the 
stiffness  of  the  arch  decreases. 
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Figure  C-8.  Load-time  relationship. 


The  cross-sectional  properties  of  the  arch  in  the  elastic  range  are  represented  by  the 
stiffness  parameters,  EA  and  El,  as  follows: 


EA  *  177.4xl06  lb 
El  =  407. 1  x  106  lb-in.2 


The  cross-sectional  properties  of  the  arch  in  the  inelastic  range  are  more  complicated; 
the  actual  cross  section  is  represented  by  two  flanges.  Attempts  were  made  to  provide  the 
closest  approximations  of  the  actual  cross-sectional  properties.  The  following  set  of 
parameters  were  used: 


Area  of  top  flange  A^ 

Area  of  bottom  flange  Aj^ 

Modulus  of  elasticity  of  top  flange  E{ 

Modulus  of  elasticity  of  bottom  flange 
Yield  strength  of  top  flange  f  ^ 

Yield  strength  of  bottom  flange  f  b 

Distance  between  flanges  u 


3.815  in.2 
2. 288  in. 2 
30  x  106  psi 
28  x  106  psi 
55,  300  psi 
66,  500  psi 

o  i  e  4  — 

1U  111* 


In  the  absence  of  specific  information  on  the  damping  characteristics  of  the  arch 
material,  it  was  assumed  that  the  external  damping  measured  by  experimental  means  repre¬ 
sented  the  total  damping  of  the  system.  An  average  damping  factor,  v ,  of  0.  095  was 
observed  from  the  fundamental  antlsymmetrlcal  mode  of  the  free  vibration  tests.  In  order 
to  conform  with  the  damping  term  used  in  the  computer  program,  the  above  damping  factor 
was  converted  to  that  of  a  complete  ring  oscillating  at  the  frequency  corresponding  to  the 
"breathing"  mode. 
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As  shown  in  Appendix  D,  the  first  antisymmetrical  frequency  of  the  arch  was  59.  9  cps, 
and  the  "breathing"  mode  frequency  of  a  complete  ring  was  317.3  cps.  The  equivalent 
damping  factor,  v',  was  as  follows: 


v' 


(0.095)  =  0.018 


Presentation  of  Numerical  Results 

Dynamic  responses  of  the  arch  in  the  elastic  range  were  obtained  for  peak  loads  of 
200  lb  'in.  and  400  lb  'in.  For  loads  above  400  lb  in. ,  the  inelastic  behavior  of  the  arch  was 
anticipated  and  the  inelastic  program  was  used.  As  a  comparison,  the  load  of  400  lb  in. 
was  repeated  in  the  inelastic  solution.  The  loading  was  increased  in  increments  of  100  lb  in. 
up  to  700  Ib/in.  and  then  increased  in  50-lb  in.  increments.  A  maximum  dynamic  load  of 
800  lb  in.  was  used  in  the  solutions. 

The  results  of  dynamic  solutions  for  elastic  response  for  peak  loads  of  200  Ib/in.  and 
400  lb/in.  are  given  in  Table  C-l,  along  with  the  maximum  values  of  bending  moments, 
axial  forces,  radial  and  tangential  displacements,  and  the  corresponding  times  when  the 
maximum  values  occurred.  The  inelastic  responses  for  loading  between  400  lb  in.  and 
800  lb/ln.  are  given  in  Tables  C-2,  C-3,  and  C-4.  AU  solutions  were  based  on  a  16-bar 
system;  however,  data  are  presented  for  intervals  of  one-eighth  of  the  arch  length.  At  the 
maximum  load  of  800  lb  in. ,  the  convergence  of  radial  and  tangential  displacements  was 
slow.  The  maximum  displacements  were  not  obtained  at  several  locations  of  the  arch  as 
indicated  in  Table  C-4. 

The  maximum  elastic  and  inelastic  values  of  moments,  axial  forces,  and  radial  and 
tangential  displacements  at  the  quarter  points  and  the  crown  of  the  arch  are  presented  in 
graphic  form  in  Figures  C-9  through  C-12.  The  bending  moment  curves  for  the  quarter 
points  coincide  with  each  other  in  the  elastic  range,  but  their  signs  are  different.  Bending 
moments  at  the  crown  are  comparative  I.,  low.  The  maximum  values  of  axial  forces  do  not 
vary  appreciably  at  the  quarter  points  and  at  the  crown  of  the  arch,  as  shown  in  Figure  C-10. 
The  maximum  radial  displacements  occur  at  the  3/4-point  and  the  maximum  tangential 
displacements  occur  at  the  crown  of  the  arch. 

The  time  interval  within  which  the  maximum  values  of  various  response  functions  are 
attained  is  only  a  very  small  fraction  of  the  entire  load  duration.  The  decrease  of  the  load 
from  its  peak  value  is  very  small  within  this  time  interval. 

The  time-variant  traces  of  moment  and  displacements  show  that  the  major  contribution 
came  from  the  first  antisymmetrical  mode  of  vibration.  This  observation  can  be  made  from 
the  moment  curves  and  displacement  curves  of  Figures  C-13  and  C-14  for  a  peak  load  of 
200  lb.  in.  The  period  of  major  oscillations  corresponds  to  the  period  erf  the  first  antisymmet¬ 
rical  mode  of  the  arch  in  free  vibration.  It  should  be  noted  that  the  influence  of  the  third 
symmetrical  mode  is  also  present,  but  the  effects  are  comparatively  small. 

The  predominant  oscillation  of  the  axial  forces  occurs  at  a  frequency  corresponding  to 
the  second  symmetrical  mode  ol  vibration,  which  is  close  to  the  "breathing"  mode  for  a 
complete  ring.  The  time  variant  traces  of  axial  forces  near  the  crown  and  supports  of  the 
arch  are  shown  in  Figure  C-15.  It  can  be  seen  that  the  variation  of  axial  forces  along  the 
arch  is  small. 

In  the  inelastic  range,  the  response  frequencies  and  modes  of  various  functions  are 
modified.  Increases  in  the  natural  period  of  the  response  curves  are  expected  as  the 
stiffness  of  the  arch  decreases. 


Table  C-l.  Maximum  Response  From  Dynamic  Loading  (Elastic  Program) 


Location 

Moment 
(in.  -lb) 

Time 

(msec) 

Axial 

Force 

(lb) 

Time 

(msec) 

Radial 
Displace¬ 
ment 
(in. ) 

Time 

(msec) 

Tangential 
Displace¬ 
ment 
(in. ) 

Time 

(msec) 

po  = 

200  lb/in. 

Left  support 

0 

0 

o 

1/8  point 

-124,000 

9.0 

-19,100 

1.5 

0.202 

8.9 

8.3 

1/4  point 

-161,000 

8.5 

-19,500 

8.0 

0. 278 

8.6 

-0. 080 

8.3 

3/8  point 

-116,000 

7.9 

-19,100 

8.0 

0.187 

7.4 

-0. 134 

8.2 

Crown 

30,  000 

3.2 

-18,800 

1.5 

10.3 

-0.155 

8. 1 

5/8  point 

118,000 

9.0 

-18,700 

1.5 

-0.  221 

8.0 

-0. 131 

8.0 

3/4  point 

161,000 

7.2 

-19, 100 

1.7 

-0. 302 

7.8 

-0.  075 

wan 

7/8  point 

123, 000 

-19,800 

1.8 

-0.215 

7.8 

-0.  021 

mm 

Right  support 

0 

m 

1.8 

0 

0 

HI 

PQ  =  400  lb/in. 

Left  support 

0 

0 

0 

1/8  point 

-254, 000 

9.0 

-38,100 

1.4 

f  >1Ti 

8.9 

8.3 

1/4  point 

-328,000 

8.5 

-39,200 

8.0 

8.6 

-0. 165 

8.4 

3/8  point 

-233, 000 

8.0 

-38,400 

8.0 

7.4 

-0.  272 

8.3 

Crown 

61,000 

3.2 

-37, 500 

1.5 

10.2 

-0.314 

8.2 

5/8  point 

242, 000 

9.0 

-37,400 

1.5 

SBhjJK  Wk 

8.0 

-0. 368 

8.1 

3/4  point 

327,000 

-38,100 

1.7 

-0.617 

7.9 

-0.151 

8.0 

7/8  point 

250,  000 

Kftl 

-39,600 

1.8 

-0.438 

7.8 

-0.041 

7.1 

Right  support 

0 

H 

lgTi]ES!i» 

1.8 

0 

0 

Table  C-2.  Maximum  Response  From  Dynamic  Loading  (Inelastic  Program) 


Location 

Moment 
(in.  -lb) 

_ 

Time 

(msec) 

Axial 

Force 

(lb) 

■ 

Radial 

Displace¬ 

ment 

(in.) 

Time 

(msec) 

Tangential 

Displace¬ 

ment 

(in.) 

Time 

(msec) 

P0  =  400  lb/in. 

Left  support 

0 

0 

0 

1/8  point 

-254,000 

9.1 

^KTjWjTjTS 

1.5 

9.0 

SEX3&B: 

8.4 

i/4  point 

-326.000 

8.5 

8. 1 

0.  552 

8.6 

-0. 162 

8.4 

3/8  point 

-234,000 

8.0 

-38,400 

8. 1 

7.2 

8.3 

Crown 

61,000 

3.3 

-37,500 

1.5 

-0.  062 

10.2 

-0.312 

8.2 

5/8  point 

242,000 

9.1 

-37,400 

1.5 

-0.454 

8. 1 

-0. 263 

8.0 

3/4  point 

326,000 

7.2 

-38, 100 

1.8 

-0.614 

7.9 

-0.151 

mom 

7/8  point 

251,000 

7.8 

-39,700 

1.8 

-0. 434 

7.8 

-0.041 

Esa 

Right  support 

0 

0 

_ 

0 

mm 

P0  =  500  lb/in. 

Left  support 

0 

mi 

0 

1/8  point 

-287,000 

9.2 

-47,500 

1.5 

■ 

9.2 

-0.  066 

8.5 

1/4  point 

-371,000 

6.4 

-47,000 

1.5 

9. 1 

8.7 

3/8  point 

-269,000 

6.6 

-47, 000 

1.5 

0.461 

7.0 

-0.350 

8.6 

Crown 

95,000 

10.2 

-46,800 

1.5 

-0. 122 

10.5 

8.5 

5/8  point 

307,000 

9.2 

-46,700 

1.5 

-0.  580 

8.4 

8.2 

3/4  point 

408,000 

7.8 

-47,600 

1.8 

-0. 783 

7.9 

-0. 190 

7.5 

7/8  point 

306, 000 

■M 

-49,500 

1.8 

-0.  546 

7.9 

-0.050 

7.2 

Right  support 

0 

wm 

-50,700 

1.6 

0 

0 

Table  C-3.  Maximum  Response  From  Dynamic  Loading  (Inelastic  Program) 


Location 

Moment 
(in.  -lb) 

Time 

(msec) 

Axial 

Force 

(lb) 

Time 

(msec) 

Radial 
Displace¬ 
ment 
(in. ) 

Time 

(msec) 

Tangential 
Displace¬ 
ment 
(in. ) 

Time 

(msec) 

po  = 

600  lb/in. 

Left  support 

0 

■■ 

0 

0 

1/8  point 

-301,  01k, 

9.8 

-56,900 

■39 

0.690 

9.8 

10.5 

1/4  point 

-393, 000 

6.4 

-56,400 

Cfl 

0. 991 

9.6 

-0. 296 

10.2 

3/8  point 

-278,000 

5.5 

-56,400 

EH 

0.  555 

9. 1 

9.9 

Crown 

129,000 

9.8 

-56,200 

■H 

-0.  191 

10.9 

-0. 530 

9.8 

5/8  point 

354,000 

10.4 

mm 

-0.  785 

10.6 

-0. 437 

9.  1 

3/4  point 

442, 000 

8.0 

-57, 100 

1.8 

-1.007 

8.3 

-0.  243 

8.9 

7/8  point 

330, 000 

6.0 

-59,400 

1.8 

-0.  672 

8.3 

-0.062 

9.0 

Right  support 

0 

1.8 

0 

0 

Pc  = 

700  lb/in. 

Left  support 

0 

■  m  B 

■MB 

0 

1/8  point 

-313,000 

12.9 

-66,400 

1.5 

1.010 

BUS 

-0. 137 

13.2 

1/4  point 

-408,000 

9.8 

-65,700 

:SSfl 

1.428 

eh 

-0.4&0 

13.1 

3/8  point 

-278,000 

5.1 

-65,700 

BSI 

0.756 

15.3 

-0.719 

13.  1 

Crown 

139,000 

9.8 

-65,500 

In 

-0.  261 

11.5 

-0. 792 

13.3 

5/8  point 

388, 000 

10.7 

-65,400 

Bn 

-1.  147 

14.1 

-0.656 

13.4 

3/4  point 

464,000 

8.0 

-66,500 

1.8 

-1.  524 

14.2 

-0  352 

13.3 

7/8  point 

352, 000 

5.6 

-69,300 

1.8 

-0.  960 

14.0 

-0.081 

12.9 

Right  support 

0 

-70, 900 

1.8 

0 

O 

Table  C-4.  Maximum  Response  From  Dynamic  Loading  (Inelastic  Program) 


Location 

Moment 
(in.  -lb) 

Time 

(msec) 

Axial 

Force 

(lb) 

Time 

(msec) 

Radial 
Displace¬ 
ment 
(in. ) 

Time 

(msec) 

Tangential 
Displace¬ 
ment 
(in. ) 

Time 

(msec) 

Po  = 

750  lb/in. 

Left  support 

0 

mm 

0 

o  ! 

1/8  point 

-317,000 

12.9 

-71,100 

bh 

1.314 

19. 1 

Ml  jifl 

1/4  point 

-413,000 

9.7 

-70,400 

BSI 

1.880 

19.1 

3/8  point 

-279,  000 

5.0 

-70,400 

R3B 

0.  990 

19.0 

19.3 

Crown 

148, 000 

9.7 

-70,200 

Bn 

-0.  346 

20.0 

-1.068 

19.2 

5/8  point 

399, 000 

11. 1 

-70, 000 

Bn 

-1.561 

18.0 

-0.887 

19. 1 

3/4  point 

473, 000 

14.3 

-71,300 

1.8 

-2.055 

17.8 

-0.462 

19.2 

7/8  point 

355, 000 

5.6 

-74,200 

1.8 

-1.262 

17.6 

-0.  098 

19.  Z 

Right  support 

0 

-75,900 

1.8 

0 

0 

P0  - 

800 lh/i 

nt 

Left  support 

0 

mm 

BN 

0 

o 

1/8  point 

-326, 000 

ESI 

-75,700 

iff! 

1.876 

38.8 

Mil.  1  fli 

40.0* 

1/4  point 

-426,000 

Warn 

-75,000 

RIB 

2.645 

39.0 

40.0* 

3/8  point 

- 279, 000 

RSI 

-74,900 

Bn 

1.344 

38.8 

40.0* 

Crown 

152,000 

■MM 

-74,800 

1.5 

-0.  558 

40.0* 

-1.567 

39.2 

5/8  point 

407,000 

EQEfl 

-74,600 

1.5 

-2.345 

40.0* 

-1.320 

40.0* 

3/4  point 

492,000 

BSI 

-76,  000 

1.8 

-3.062 

40.0* 

-0.657 

39.0 

7/8  point 

364,000 

BH 

-79,  200 

1.8 

-1.845 

40.0* 

-0.120 

38.8 

Right  support 

0 

HI 

-81,000 

1.8 

0 

0 

♦End  of  solution;  maximum  values  were  not  obtained. 
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(lb  X  l(f3) 


Tin 


axial  forces. 


Load  (Ib/in 


Radial  Displacement  (in.) 


Figure  C-ll.  Maximum  radial  displacement  for  dynamic  loads. 
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Figure  C-13.  Bending  moments  versus  time. 


Appendix  D 

ANALYSIS  FOR  NATURAL  FREQUENCIES  AND  BUCKLING  LOADS 


r 


In  the  analysis  of  the  two-hinged  circular  arch  for  its  natural  frequencies  and  buckling 
l  -  the  active  continuous  structure  of  Figure  B-l  is  reduced  to  a  discrete  framework 
consisting  of  a  series  of  rigid  bars  and  flexible  joints  as  described  in  Appendix  C.  The 
mass  may  either  be  considered  lumped  at  the  joints  or  distributed  along  the  bars.  A  uniform 
static  normal  pressure  may  be  introduced  in  the  frequency  equations;  i.e. ,  the  frequencies 
of  the  arch  can  be  determined  for  free  vibration  or  forced  vibration.  The  theoretical  devel¬ 
opment  of  equations  was  performed  by  Eppink  (Reference  15).  Additional  remarks  are  made 
primarily  for  clarity  and  continuity  of  this  report. 

Derivation  of  Equations 

In  considering  the  uniform  static  pressure  acting  on  the  arch,  the  following  additional 
assumptions  are  made: 

1.  The  arch  maintains  its  circular  shape  under  the  external  normal  pressure. 

2.  The  arch  undergoes  only  a  uniform  compression  under  the  external  normal 
pressure. 

3.  All  measurements  are  taken  from  the  compressed  position  of  the  arch. 

4.  The  oscillation  of  the  arch  is  small. 

5.  The  effect  of  damping  is  negligible. 

The  above  assumptions  serve  to  linearize  the  vibration  of  the  arch  under  uniform 
normal  pressure.  The  axial  force,  Pjt,  can  be  expressed  in  the  form 


Pjt  =  -|wLcot|  *  N.  (D-l) 


where  the  first  term  represents  the  force  due  to  uniform  compression  and  the  second  term 
represents  the  force  due  to  displacements.  By  substituting  Equations  C-3  and  D-l  into 
Equations  C-l  and  C-2  and  omitting  the  product  terms  of  N^  and  M^,  the  following  expres¬ 
sions  for  Zjr  and  Zjt  are  obtained: 


z)r  *  -(Ni  *  Nk)51"!  ’  CK  '  2M)  *  Mk)“s  I 

+  wLcos^  +  "2wL(^k  "  ^j)cot4cosT  (D-2) 

S  =  -(Nj  -  Nk)cosf  +  rK  -  Mk)sln  I 

-r  |wL^k  +  0.)  cos  (D-3) 
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If  the  axial  deformation  is  taken  into  account,  the  components  Wjr  and  Wjt  of  the 
external  forces  are  established  by  noting  that 

p,1  -  4(l  +  dj)w 


p 


i 


r 


w 


Substituting  the  above  expression  into  Equations  C-4  and  C-5,  and  neglecting  the  product 
terms  of  dtp,  the  following  expressions  for  Wjr  and  Wjt  are  obtained 


-wL  cos  ^ 

■  iw(*k +  K 

)cos|  + 

(D-4) 

wk  ■  4* 

K  -  dj)sinl 

-  +  ^j)cosl 

(D-5) 

The  components  of  the  internal  force,  Zjr  and  Zu,  of  Equations  D-2  and  D-3  can  be 
expressed  in  terms  of  deformations  i  and  tp  by  using  Equations  C-16  and  C-17.  The 
resulting  equations  along  with  Equations  D-4  and  D-5  are  then  expressed  in  terms  of  displace¬ 
ments  by  using  the  first-order  approximations  of  Equations  C-8,  C-9,  and  C-12.  The 
resulting  expressions  for  Zjr,  Zjt,  Wjr,  and  Wjt  are  then  expressed  in  terms  of  displace¬ 
ments  djr  and  d;t. 

At  a  natural  mode  of  vibration  the  conditions  of  simple  harmonic  motion  are  satisfied 
for  each  particle  of  the  system,  and  the  equations  of  motion,  Equations  C-6  and  C-7,  become 


j  jr 


ii.  u.. 

J  Jt 


'J  V  =  Zjr  +  Wjr 

(D-6) 

1.  =  Z  ,  +  W.„ 

1  Jt  Jt  jt 

(D-7) 

where  u  is  the  circular  frequency  of  vibration. 

Substituting  the  expressions  for  Zjr,  Zjt,  Wjr,  and  Wjt  into  Equations  D-6  and  D-7, 
simplifying  by  letting  p  denote  the  mass  per  unit  length  of  the  arch  and  r  the  radius  of 
gyration  of  the  arch  cross  section,  and  introducing  the  following  notations 


c 


(D-8) 


P 


(D-9) 
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=  (4csinl)  ^w2 


(D-10) 


k  =  4c  sin'1 


(Dtf 


(D-ll) 


the  following  expression  for  the  equation  of  motion  in  the  radial  direction  is  obtained: 

[v  +  (-4  +  pc2)\r  +  (6  +  2^c2)4ir  *  (-4  +  *c2Kr  *  4Ir] 

+  c  [4ht  -  <2  +  p)  4it  +  (2  *  p)  4kt  *  Alt] 

-  ^jr  *  k  [-(2  *  7)  4ir  *  72  V  '  (2  *  \)  \r  *  c  Hit  *  \t)]  ^ 


The  corresponding  equation  for  the  tangential  direction  becomes 
c  [~Ahr  +  (2  +  p)  4ir  '  (2  +  p)  4kr  +  4ir] 


[-c2dht  -  pd.t  .  2  (p  +  c>jt  -  pdkt  -  c2  dltj 


=  Xd 


«  *  k  [c  K 


- 

ir  kr 


)  *  K  • 2  v 


)] 


(D-13) 


The  left  sides  of  Equations  D-12  and  D-13  are  applicable  to  joints  2  <  j  <  n-2.  The 
boundary  conditions  must  be  considered  in  deriving  the  equations  of  motion  for  joints  1  and 
n-  1.  Figure  D-l  is  a  schematic  representation  of  the  equations  of  motion  for  the  general 
interior  joint  j  and  the  joints  near  the  supports. 

A  system  of  2  (n-  1)  homogeneous  linear  algebraic  equations  can  be  obtained  from 
Equations  D-12  and  D-13  with  specialized  equations  for  the  joints  near  the  supports.  The 
system  of  equations  can  be  written  in  abbreviated  form  as  lollows: 


Ad 


\  T>  A 
aDu 


r  n  A 


0 


(D-14) 


where  A,  B,  and  C  are  symmetric  matrices  of  order  2(n-  1),  and  d  is  a  ulumn  matrix 
of  the  displacement  dr  and  dt.  There  are  2  (n  -  1 )  characteristic  roots  and  vectors  of  the 
above  equations  corresponding  to  the  natural  frequencies  and  mode  shapes  of  vibration  of 
the  system. 
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Figure  D-l.  Left  side  of  Equations  D-12  and  D-13  in  block  form. 


0 


It  should  be  noted  that  for  free  vibration  the  load  parameter  k  is  set  to  zero  and 
Equation  D-14  becomes 


A  d  -  ABd  =  0 


For  the  buckling  loads  of  the  system,  the  frequency  parameter  A  is  set  to  aero  and 
the  following  expression  is  obtained: 


Ad  -  kCd  =  0 


The  mass  of  the  system  can  be  considered  to  be  distributed  along  the  bars  of  the  discrete 
framework  instead  of  concentrated  at  the  joints.  The  equations  can  be  established  by  using 
the  Lagrange  equations  of  motion.  A  detailed  derivation  is  presented  in  Reference  16.  The 
resulting  equations  can  be  presented  in  abbreviated  form  as  follows: 


Ad  -  AB'd  -  kCd  =  0 


where  B*  is  a  symmetric  matrix  of  order  2  (n  - 1).  The  terms  AB’d  can  be  obtained  for  the 
equations  of  motion  in  the  radial  direction  as  follows: 

A^-(cos  0dir  +  4 djr  +  cos  #4kr  +  sin  tfdlt  -  sin  (D-15) 

and  for  equations  of  motion  in  the  tangential  direction  as  follows: 

A^-sin  #dir  +  sin^dkr  +  cos  ^d^  -r  4d-t  *  cos  #dktj  (D-16) 


For  geometrically  symmetrical  structures,  the  natural  modes  of  vibration  are  either 
symmetrical  or  anti  symmetrical  with  respect  to  the  line  of  symmetry  of  the  structure.  This 
can  be  utilized  to  determine  the  complete  set  of  natural  modes  and  the  associated  frequencies 
by  solving  two  characteristic  value  problems,  each  of  order  (n-  1)  instead  of  one  of  order 
2  (n  - 1). 

Let  joint  k  be  located  at  the  line  of  symmetry  of  the  discrete  framework.  For 
antisymmetrical  modes,  the  radial  displacement,  dj^,  at  joint  k  is  zero,  and 


^lr  “  ”^jr  ’  ^it  4jt 


d  =  -d. 
mr  ir 


d  ,  =  d.. 
mt  it 
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For  symmetrical  modes,  the  tangential  displacement,  at  joint  j  is  zero,  and 


4 


4r 


jr 


It 


-4 


it 


d 

mr 


ir 


mt 


-d 


it 


Figure  D-2  presents  the  left  side  of  Equations  D-12  and  D-13  for  joints  j  and  k  in  block  form. 

As  an  example,  the  frequency  equations  for  a  12-bar  system  are  presented  in  matrix 
form  in  Tables  D-l  and  D-2  for  antisymmetrical  modes  and  for  symmetrical  modes, 
respectively.  In  order  that  the  same  matrix  equations  be  applicable  to  frequency  as  well  as 
buckling  load  calculations,  three  additional  parameters,  x,  y,  and  z,  are  introduced. 

For  free  vibration  of  concentrated  mass  systems 


x  =  X  ,  y  =  Q  ,  z  =  0 


For  free  vibration  of  distributed  mass  systems 


x 


y 


<t> 


-|-sin  <f> 


For  buckling  loads 


x 


z  -  -ck 


Computer  Program 

A  computer  program  has  been  completed  which  can  calculate  either  the  natural 
frequency  or  buckling  load  of  an  arch  of  uniform  cross  section  based  on  the  equations  given 
in  Tables  D-l  and  D-2.  The  source  program  is  written  in  FORTRAN  and  can  be  carried  out 
on  an  IBM  7094  or  IBM  16201!  computer.  It  should  be  noted  that  the  number  of  bars  consid¬ 
ered  is  limited  only  by  the  memory  capacity  of  the  computers.  However,  solutions  for  more 
than  40  bars  are  not  recommended  considering  the  computer  time  involved  and  the  very  slight 
increase  in  accuracy. 

The  natural  frequency  can  be  determined  for  either  a  concentrated  mass  system  or  a 
distributed  mass  system.  For  natural  frequency  calculations,  a  uniformly  distributed  load 
can  be  acting  on  the  arch.  A  special  case  occurs  when  the  natural  frequency  vanishes:  the 
buckling  load  for  the  system  is  obtained. 

The  source  program  can  be  subdivided  into  the  following: 

a.  Determination  of  natural  frequencies  and  buckling  load 

b.  Determination  of  natural  modes 

It  is  possible  to  divide  the  program  so  that  a  smaller  computer  can  be  used  for  the  solutions. 
To  facilitate  this  operation,  complete  flow  charts  are  given  in  Figures  D-3,  D-4,  D-5,  and 
D-6  for  general  analysis  of  arches,  evaluation  of  matrices,  search  for  roots,  and  vector 
components.  The  source  program  in  FORTRAN  is  given  in  Program  D-l. 


Antisymmetry  1  H  -4+pc^  H  5+2pe2  H  -o  *  c  h  -c(2+p)  H  -c  H  c(2+p) 
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Figure  D-2.  Left  side  of  Equations  D-12  and  D-13  in  block  form. 


Table  D-2.  Frequency  Equations  of  Symmetrical  Modes  for  Circular  Arches  in  the  Elastic  Range 


START 


Read  control  N. 
cord  w/  Nroot, 
Xo,  Fmax,  Delta, 
Tol,  Xk  J 


Solve 

homogeneous 

simultaneous 

equations 


^  Fmax- 
Abs  (F2) 


lretf-1 


-  .  . 

Raise  X2 
by  Delta 

Determinant 

A 

Search  for 
root  by  nested 
interval  method 


Compute 

residual 

vector 


Output 
solution  and 
residual  vectors 


Output 

root,  residual, 
frequency  and 
period 


Figure  D-3.  Flow  diagram,  analysis  of  arches. 


Matrix  A  “  0 


Ifreq  -5 


Matrices  for 
frequencies  and 
rrode  shapes 


Matrix  for 
buckling  load, 
antis/  mmetrica! 


Elements  of 
matrices  A, 
concentrated 
mass  system 


Elements  of 
matrices  A, 
antis/mmetrieal 
modes 
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Program  D-l.  Natural  Frequency  and  Mode  of  Circular  Arch 


901 

902 
100C 


loot 


1002 

1003 
2000 

C 

1004 


1005 

1009 

1010 
1011 

1013 


1015 


program  of  circular  ARCH  ANALYSIS 

purpose  of  program*  to  evaluate  x2  in  matrix  of  ax*o« 

program  DETERMINES  NATURAL  frequencies  of  CIRCULAR  arch 

ANTI-SYMMETRICAL  AND  SYMMETRICAL  MODE  SHAPES 

BUCKLING  LOAD  AND  DEFLECTED  SHAPE  OF  THE  ARCH 

PROGRAM  CONTRuLLED  BY  VALUES  OF  IFREQ 

JFREU  1,  CONCENTRATED  MASS  SYSTEM  ANTI-SYMMETRICAL  MODE 
IFREQ  2.  CONCFN  TR  A  f ED  MASS  SYSTEM  SYMMETRICAL  MODE 
I  F'REQ  3  .  DISTRIBUTED  MASS  a  Y  S  T  FM  ANT  I -S  YMME  TR  I  CAL  MODE 
I  FREQ  4,  DISTRIBUTED  MASS  5  Y  ST  E v  SYMMETRICAL  MODE 
IFREQ  5  ,  BUCKLING  LOAD  OF  AN T I -SYMMf TR I  CAL  MODE 

ifreq  6,  buckling  load  of  symmetrical  mode 

MOOIFIED  25  JUNE  1965  I  RICHARD  H.  CHIU  > 

DOUBLE  PRECISION  A  ,  X  «  R  «  DUMP  *DET  .FO.F1.F2.F3 
DIMENSION  A ( 12,12 ) .XI 12  > 

FORMAT  (  I5.E15.8.D15.8.3E15.8> 

FORMAT  I / 1 1 X  .4HROOT  .9X.8HRESI DUAL  «BX,5HFREQ.«10X,6HPER|0D/) 

OC  1032  I F  RED* 1*6 

READ  INPUT  TAPE  5.  CONTROL  CARD 

READ  15.  1001  INROOT .XO.FMAX, DELTA. IOL.XK 

FORMAT! 15.5E15.8l 

READ  INPUT  TAPE  5.  INPUT  PARAMETERS 

N*ORDER  OF  DETERMINANT 

RC  =  RADIUS  OF  ARCH 

PHl«  INCREMENT  OF  CENTRAL  ANGLE 

RL=  RADIUS  OF  GYRATION  OF  ARCH  SECTION 

U«MU 

EI«E«I 

READ  (6.  10011  N.RC.PHI .RL.U.EI 

XL»2.*RC*SIN  (PHI/2.) 

C2 *S I N  (PHI/2. >/CDS  TPhI/2.) 

FOR  BUCKLING  LOAD 

IF  ( I  FREQ-5 )  1003.1002,1002 

C1«(C2*XL**3)/(2.«£I ) 

GO  TO  2000 

Cl«XL**4*u/(EI*COS  (PHl/2. )**2) 

RHO»(XL/RLI**2 
DO  1032  JROOT *  1 .NROOT 
WRITE  (6.  902) 

BEGIN  SEARCH  FOR  ROOT 

I  RE  T  *  1 

I  RE  TF  *0 

X1=X0 

X2*Xl 

GO  TO  1016 
F 1  =F2 

WRITE  (6,  9u 1 )  IRET.X2.F2 

I  RE  T  *2 
X2*X1 +DEL  T  A 
GO  TC  1016 

IF(F1*F2)1013, 1013, 1011 

X1-X2 

F1=F2 

WR I i -  <5.  9u  1  )  IRtT.X2.Fz 

GO  TO  1009 
X  3  *X2 
F3  =  F2 


X4  =  X3 
I  RE  T  *  3 


XO  *X  3 
F0«F3 

WRITE  (6.  901) 

X2  =X 1  +  { X  3-X 1 )»F1/(F1-F3) 

GO  TO  1016 


1RET.X2.F2 


Continued 


116 


nAAn  n  n  a  a  aa  aa  a  aaa 


Program  D-l.  Natural  Frequency  and  Mode  of  Circular  Arch  (Contd) 


evaluation  of  determinant  of  order  n 
A.  PLOT  K I N  AND  W,  W  t  L  COXSON 

29  AUGUST  1963 

3005  0E T» 1 • 

00  3013  *«2.N 
KM1«K-1 

search  for  NON-ZERO  PIVOT  ELEMENT 
DO  3006  K1*(CM1,N 
I F  <  A ( KM  1 ♦ Kl I ) 300  7 .3006.300  f 

3006  CONTINUE 
DETERMINANT  IS  ZERO 
0ET*0. 

GO  TO  30 1 A 

DIVIDE  COLUMN  bY  PIVOTAL  tLEMEM 

3007  DO  3008  K2*K.N 

3008  AIS2.K1 I «A( K2.R1 I /AdCMl ,K1 » 

IFIKMl-Kl I30C9.3011.3C09 

MOV  t  PIV01  ELEM.  to  UPPER  LEFT  HAND  CORNER  OF  REOUCtD  DE I E  RM I  NAN  I 
UPDATE  SIGN  F  OR  INTERCHANGING  COLUMNS 

3009  DETs-DEI 
DO  3C10  X2*XMl»N 
DUMP«A(K2.KM1 I 
AIX.2.KM1  )*AIK2.K1  I 

3010  AIK2.K1)  “DUMP 

REDUCE  ORDER  OF  DETERMINANT  BY  ONE 

3011  DO  3012  I=K.N 
DO  3012  J  *  K  »  N 

3012  A<  l  .JI*A(  I  ,J)-AI  I  .KM1  )*A<r.Ml  ,J] 

3013  DET=DET»A(KM> ,KM1 ) 

3014  F2=DET*A!N.N) 

GO  TO( 1005.1010.10171 ,1RET 

SEARCH  FOR  ROOT  BY  METHOD  OF  NESTED  INTERVALS 

1017  iFIAbS  <  X*.-X2)-TOL  )  1029. 1029. 1018 

1018  IFIX3I1019. 1020. 1019 

1019  IFIAbS  I  l.-X2/X'.)-IOL)1029. 1029. 1020 

1020  IF(F1*F3)1021, 1021. 1029 

1021  IF(F1«F2)1026. 1029. 1025 

1025  X1=X2 
X4«X1 

F  1  *F2 

GO  TO  1015 

1026  X 3 *X 2 
X4*X3 
F3  =  F2 

GO  TO  1015 

IGNORE  IF  ASYMPTOTE 

1029  IFIFMAX-ABS  I F2 I > 1 004 . 1 004, 103 0 

1030  XF*X2*X2 
P=1./XF 

WRITE  16.  901 ) JROOT.X2.F2.XF  ,P 

GO  TO  lu 16  TO  EVALUATE  ELEMENTS  OF  COEFFICIENT  MATRIX 
I  RE  T  F  *  1 
GO  TC  1016 

SUb-PROGRAM  TO  SOLVE  HOMOGENEOUS  S I MuL i ANtuuS  EOuATIUNS  ISANK.-N— li 
W.  WILCOXSON 

30  august  1963 
OUTPUT  M4TRU  A 

7000  WRITE  16.  75.0) 
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7500  FORMAT  I 14X.8HMATR1X  A) 

DC  3  JJ*1,N 

DO  3  T  T  *  1  • 

3  WRIT':  (6*  2) JJ.1I .AiJj.lI) 

2  FORMAT <2HA<  .I3.1H,I3.2H)«D14.8) 

C  NO  1 6  .’HAT  N-MZ*RANK  OF  COEFFICIENT  MATRIX 

MZ=C 

00  7028  Ml*l.N 
M5«M1-MZ 

C  SEARCH  FOR  NON-ZERO  PIVOT  ELEMENT 

00  70C 1  M2*M5»N 

I F ( TOL-ABS  I  At  M2* Ml m  7002*7001  ;7001 

7001  CONTINUE 
MZ=MZ+1 

GO  TO  7028 

7002  I F ( M2-M5 )7003»7U05»7G03 

C  INTERCHANGE  ROWS  TO  OBTAIN  NON-ZERO  PIVOT  ELEMENT 

7003  DO  7004  M3=M1,N 
DUMP*A|M5,M3) 

A(M5,M3)*A(M2,M3l 

7004  AIM2.M3)=DUMP 

7005  DO  70Co  M3*M1,N 
M2*N-M3+M1 

7006  A(M5.M2I*A(M5,M2l /AIMS. Ml  I 
DO  7018  M4»1,N 

|f  (M4-M?. )  7007, 7018. 7u07 

7007  DO  7008  M3=M1,N 
M2*N-M3+M1 

7008  A(M4,M2)=A(M4,M2I-A(M4,M1 )*A!M5»M2> 

7018  CONTINUE 

7028  CONTINUE 

C  WRITE  OUTPUT  TAPE  6.  TOLERABLE  ZERO 

WRITE  (6,  5  I  A !  N .N ) 

5  FORMAT! 15HTOLERABLE  ZERO*D15, 8) 

C  VERIFY  THAT  RANK  IS  N-l 

IF (MZ-1 J 7009.7013 .7009 

7009  M1=N-MZ 

WRITE  ( 6 .  7010)  Ml 

7010  FORMAT  (  11HERR  RANK*'.  I  41 

7011  WRITE  <6.  7u 12 ) 

7012  FORMAT < 14HZER0  COMPONENT) 

C  VERtFY  THAT  ALL  SOLUTIONS  ARE  NON-ZERO 

7013  NM 1 =N— 1 

DO  7014  M 1  *  1  .  NM1 

Ir < 1.-ACM1, Ml ) »701 1,7014,7011 

7014  X  I  M 1  )  =  —  A I M 1  *  N ) 

X  C  N )  =  1 » 

CALCULATE  RESIDUALS  BY  MATRIX  MULTIPLICATION 
EVALUATE  ELEMENTS  OF  MATRIX 
I  RET  F  =  2 
GO  TO  1016 

5004  WRITt  (6,  75ul) 

7501  FORMAT!  17H  SOLUTION  VLC T OK • 1 5X . 1 5HR t S I  DUAL  VtCTOR! 
DO  5006  I  1 1 *N 

R*0, 

DO  5005  J* 1 ,N 

5005  R=AI I,JI*X! J)+R 

C  OUT  PjT 

5006  WRITE  (6.  5uG7  )  I  .  X I  I  )  ,  I ,R 

5007  FORMAT ( 2HXI  ,  1 2 ,2H )  =  , D 14.8 « 1  OX « 2HR (  , I i , 2H )  =  * 014,8 ) 
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lu32  COM  i sue 

GO  TO  10''0 

c  USE1?  MUST  I  NsERT  FUNCTION  EVALUATION  STATE  MINTS  AFTER  1016 

>:  AND  THEN  COMPILE  AN  OBJtCT  PROGRAM 

C  EVALUATION  ROUTINE  FUR  MATRIX  ELFMENTo 

1016  CONTINUE 

00  2017  1*1  .N 
00  2017  J*I  .N 

2017  A( l »J) >0. 

IF  (  I  FREQ-5 )  2ul9.2ol8.20lb 

2018  ALPHA«C1*X2 
GO  TO  2020 

2019  ALPHA*C1*39.478417*X2**4 

2020  N1*N-1 
NW*Nl/2 

All  ,1 )*5.+2.*RHO*C2*C2 
J 1 *NW+ 1 
A ( 1 • J1 ) *C2 

A ( J 1  ,  J1 ) =2. *RHO+C2*C2 
GO  TO  (  121.122.121.122.121 .1221  ..FREQ 
C  ANTI-SYMMETRICAL  MOOES 

121  A(NW.NW)=A( 1.1) 

AlNM.Ni ) *-C2 
A(N1 ,N1 I = A ( J 1 . J1 ) 

A ( NW-1 ,N)*-C2 
A(NW.Nj=C2*( 2 .+RH0 ) 

A(N-2»N)*-C2*C2 

A ( N1 ,N ) *-RH0 
A(N.N).A<J1.J1 J*C2*C2 
GO  TO  123 

C  SYMMETRICAL  MODES 

122  A ( NW  *Nw I  - A ( ! .1 1  +  2. 

A(NW.Nl)*  C2 

A[Nl,Nl)«A<Jl.Jll+2.+C2*C2 
A ( NW— 1 .N ) s 1 . 

AtNW.N) «RH0*C2*C2-4. 

AIN-2.N1*  C  2 

A(N1»N1«-C2*(2.+RH01 

4<N»N)*All«ll+l. 

123  1=1 

00  101  J*2,NW 
A  I  l  ,J)»RH0*C2*C2-4. 

101  I  *  I  ■*- 1 
1*1 

00  102  J *3 »NW 
A  (  I  .  J  )  *  1  * 

102  1*1+1 
J2  *Nw- 1 

00  106  J=2,J2 
105  A ( J  ,  J ) =  A ( 1 , 1  |  +  1. 

1*1 

J1 *NW+2 

DO  103  J 1 *N1 
A (  I.J)=C2*(2.+Rh0] 

103  1*1+1 
1*1 

Jl*NW+3 

00  104  J*J1.N1 
A  (  I  .  J ) *-C2 
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104  I-I  +  l 
1-2 

Jl-NW+l 

J2-N-2 

00  106  J-J1.J2 
A  (  I  *J)--C2*(2.-+RHO) 

106  1=1+1 
1-3 

J2-N-3 

00  107  J-J1.J2 
A(  l  .  J  )  -C2 

107  l-l+l 
l-NW'+l 
Jl-NW+3 

00  108  J-J1.N1 
A <  I.J)  —  C2*C2 

108  I-I  +  l 

I  -NW-+1 
J1 -NW+2 

DO  109  J-Jl.Nl 
A  (  1  • J) — RHO 

109  I-l  +  l 
Jl-NW+2 
J2-N-2 

DO  110  J-J1.J2 

110  A(J.J!«2.*(RHO+C2*C2i 
NW-NW+1 

GO  TO  <  125. 125. 131. 131. 140.140), IFREQ 
C  BUCKLING  LOAD 

140  ASN=-ALPHA»C2 
DO  143  J-NW  .N 

143  A< J. JJ«A( J.J>-2. -ALPHA 
J1-NW-+1 

DO  144  J-J1.N 

144  A<  J-1.J)-A<  J-1.JI-+ALPHA 
ACS— ALPHA*(2.  +  1./(C2*C2|  I 
ALPHA =ALPHA*2»/(C2*C2 ) 

J2-Nta-1 

DO  141  J-2.J2 

141  A< J-l.J) =A( J-l ,J)-ACS 
DO  142  J-1.J2 

142  A) J,J)=A( J. J)-ALPHA 
GO  TO  139 

C  CONCENTRATED  MASS  SYSTEM 

125  DO  126  J-l.N 

126  A(J«J)-A(J«J) —ALPHA 
GO  TO  127 

C  DISTRIBUTED  MASS  SYSTEM 

131  ASN-ALPHA-SIN  (PHI)/6. 

ACS*ALPHA»COS  (PHI) fb. 
ALPHA=ALPHA*2./3. 

DO  132  J-l.N 

132  A(J»J)«A(J.J)— ALPHA 
J2-NW-1 

DO  133  J-2.J2 

133  A ( J-l ,J) -A( J-l ,J) -ACS 
J1-NW-+1 

DO  134  J-J1.N 

134  A(J-1*J)»A(J-1.J)-ACS 
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139  J2*N-2 
1*2 

00  133  J*NW» J2 
A ( I »  J ) * A ( I « J  )  -ASN 

135  1*1+1 
J1*NW+1 
1*1 

DO  136  J«J1.N 
A  t  I  «  J  )  *  A.  (  I  *  J ) + ASN 

136  1*1+1 

127  00  111  l«2tNl 
NC* 1-1 
00111J*1 .NC 

111  A« 1 ,JI*AC J,1  ) 

DO  112  J*1,N1 

112  A(N,J»*2.*A« JtNJ 

IF  ( IRETF-1 )  3005»7000,5C04 
END 
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Presentation  of  Numerical  Results 

The  input  parameters  used  in  the  numerical  analysis  of  the  natural  modes  and  the 
associated  frequencies  of  the  two-hinged  circular  arch  are  the  same  as  those  used  in  the 
elastic  program  of  Appendix  C.  In  the  solutions  presented  in  this  section,  the  effects  of  the 
number  of  bars  and  the  assumptions  regarding  the  mass  distribution  of  the  discrete  frame¬ 
work  are  studied. 

Solutions  for  increasing  the  number  of  bars  from  8  to  40  in  the  discrete  framework, 
considering  the  mass  to  be  concentrated  at  the  joints,  are  tabulated  in  Table  D-3.  The 
corresponding  results  for  the  distributed  mass  system  are  given  in  Table  D-4.  The  first 
six  natural  frequencies  of  the  antisymmetrical  and  symmetrical  modes  are  presented  in  both 
tables. 

A  comparison  of  the  antisymmetrical  modal  frequencies  for  the  concentrated  mass  and 
distributed  mass  systems  is  presented  in  Figure  D-7.  The  frequencies  were  plotted  against 
the  number  of  bars  used  in  the  discrete  framework.  A  corresponding  comparison  of  the 
symmetrical  modal  frequencies  is  presented  in  Figure  D-8.  It  should  be  noted  that  both 
assumptions  regarding  the  mass  distribution  of  tne  framework  lead  to  the  same  asymptotic 
frequency.  However,  the  asymptotic  frequency  can  be  obtained  much  more  rapidly  by  the 
distributed  mass  assumption. 

As  shown  in  Figures  D-7  and  D-8,  nearly  all  the  asymptotic  frequencies  are  obtained 
for  fewer  than  40  bars.  For  practical  purposes,  the  solutions  of  the  40-bar  distributed  mass 
system  may  be  regarded  as  "exact"  solutions  for  the  arch  considered  in  this  report. 

The  normalized  natural  modes  of  vibration  of  a  40-bar  distributed  mass  system  are 
given  in  Figures  D-9  and  D-10  for  antisymmetrical  and  symmetrical  modes,  respectively. 

It  is  observed  that  the  normalized  modes  are  predominantly  flexural.  However,  the  axial 
deformation  of  the  second  antisymmetrical  mode  is  significant  and  is  closely  related  to  the 
"breathing"  mode  of  a  complete  ring  having  the  same  cross-sectional  properties.  A 
comparison  can  readily  be  made  by  considering  the  frequency  of  vibration  of  a  complete  ring  of 
uniform  cross  section  as  given  by  the  following  expression: 


which  is  calculated  to  be  317. 3  cps  for  the  cross-sectional  properties  of  the  arch.  This  value 
is  comparable  to  the  frequency  of  301. 1  cps  for  the  second  symmetrical  mode  of  the  arch. 
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Table  D-3.  Natural  Frequencies  and  Periods  of  a  Two-Hinged  Arch 

Concentrated  Mass 


Number  of  Bars 


Frequencies  for  Antisymmetrical  Modes  (cps) 


Period 

T* 


Frequencies  for  Symmetrical  Modes  (cps) 


135.2 

137.5 

138.6 

139.2 

139.8 

140.1 

295.8 

298.1 

299.3 

299.7 

300.2 

300.4 

388.0 

411.5 

422.8 

429.3 

435.5 

438.8 

651.9 

739.2 

782.6 

806.9 

831.6 

843.1 

886.5 

1,099 

1,191 

1,206 

1,211 

1,213 

1,197 

1,208 

1,233 

1,286 

1,350 

1,382 

128.5 

286.1 

329.2 

414.3 

1,172 


Natural  period  (msec)  based  on  40-bar  solutions 


7. 14 
3.33 
2.  28 
1.19 
0.824 
0.  724 


Table  D-4.  Natural  Frequencies  and  Periods  of  a  Two-Hinged  Arch 


Distributed  Mass 


Number  of  Bars 

Period 

T* 

Mode 

8 

12 

16 

20 

24 

32 

40 

Frequencies  for  Antisymmetrical  Modes  (cps) 

1 

59.83 

59.86 

59.  86 

59.86 

59.86 

59.86 

59.86 

16.71 

V 

2 

270.2 

272.0 

272.3 

272.5 

272.5 

272.5 

272.5 

3.63 

V 

h 

3 

583.1 

607.1 

608.5 

608.7 

608.7 

608.7 

608.7 

1.64 

*=• 

4 

677.5 

685.4 

685.6 

685.4 

685.2 

684.7 

684.7 

1.46 

- 

5 

— 

1,095 

1,125 

1,131 

1,132 

1,133 

1,134 

0. 882 

6 

1,569 

1,728 

1, 752 

1,757 

1,759 

1,759 

0. 569 

Mode 

Frequencies  for  Symmetrical  Modes 

(cps) 

H 

140.3 

140.5 

140.5 

140.5 

140.5 

140.6 

140.6 

7.11 

i; 

ci 

302.3 

301.6 

301.4 

301.1 

301.1 

301. 1 

301.1 

3.32 

< 

B 

432.1 

442.2 

443.4 

443.6 

443.8 

443.8 

443.8 

2. 25 

B 

733.0 

849.6 

860.4 

862.8 

863.4 

863.9 

864.1 

1.16 

f- 

T 

s 

--- 

1,233 

1,226 

1,222 

1,220 

1,218 

1,217 

0.822 

6 

1,257 

1,347 

1,420 

1,433 

1,436 

1,438 

1,439 

0.695 

♦Natural  period  (msec)  based  on  40-bar  solutions. 
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Symmetrical  Modal  Frequencies  (cps) 


first  anti  symmetrica  I  mode  of  vibration,  fn  ■  60  eps 


third  antisymmetricol  mode  of  vibration,  fn  *  609  cps 


Figure  D-9.  Antisymmetrical  modes  of  two-hinged  arch. 
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Appendix  E 

SIMPLIFIED  ANALYSIS  OF  A  CIRCULAR  ARCH 


In  the  static  analysis  (Appendix  B)  of  the  two-hinged  arch  of  Figure  B-I,  the  distributions 
of  moments  and  stresses  are  basically  antisymmetrical  and  the  axial  force  along  the  arch 
does  not  vary  appreciably.  In  the  dynamic  analysis  of  Appendix  C,  the  response  of  the  arch 
is  largely  contributed  by  the  first  antisymmetric  mode  of  vibration;  the  axial  force  also 
remains  essentially  uniform  along  the  arch.  The  findings  suggest  that  a  simplified  analysis 
and  design  procedure  similar  to  those  proposed  by  Newmark  can  be  used.  All  notations 
in  this  appendix  conform  to  the  previous  definitions  unless  otherwise  noted. 

Static  Analysis 

Assuming  that  the  principle  of  superposition  holds  true  in  the  elastic  range  of  the  arch 
response,  the  uniform  load  distributed  over  one -half  the  arch  length  can  be  obtained  by  two 
components  as  follows: 

1.  Symmetrical  loading  uniformly  distributed  over  the  entire  arch  span  with  a 
magnitude  equal  to  one-half  the  original  load. 

2.  Antisymmetrical  loading  consisting  of  an  inward  uniform  loading  distributed  over 
one-half  the  arch  length  and  an  outward  uniform  loading  over  the  other  half  of  the 
arch  length.  The  magnitude  of  the  antisymmetrical  loading  is  also  one-half  the 
original  load. 

The  sum  of  the  two  components  is  equivalent  to  the  original  loading  condition. 

For  a  circular  arch  with  uniform  cross  section  along  its  length,  the  critical  section  is 
at  the  quarter  point  of  the  loaded  half  of  the  arch.  The  static  axial  force  at  the  critical 
section  due  either  to  symmetrical  or  antisymmetrical  loading  is  given  by  the  approximate 
expression  as  follows: 


Pt  =  wR 


(E-l) 


For  the  symmetrical  loading,  the  maximum  bending  moments  at  the  quarter  points  are 
given  by  the  approximate  expression 


M 


1 


(E-2) 


For  the  antisymmetrical  loading,  the  maximum  bending  momeut  at  the  quarter  point 
under  the  inward  loading  is  given  by  the  approximate  expression 


» 

»  . 
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and  the  corresponding  buckling  load,  wcr, 


can  be  obtained  from 


3 

»■ 


(E-i) 


Stress  amplification  as  a  result  of  the  tendency  of  the  arch  to  buckle  under  the  axial  force  can 
be  expressed  by  the  amplification  factor 


A 


m 


(E-5) 


The  amplification  factor  will  be  used  as  a  multiplier  for  the  moment  terms  of  Equations  E-2 
and  E-3.  The  total  static  moment  is  the  sum  of  Mi  and  M2;  the  stresses  and  strains 
caused  by  the  bending  moment  and  the  axial  force  can  be  obtained  by  applying  Equations  B-30 
through  B-33.  The  static  yield  condition  may  be  established  with  the  aid  of  stress-strain 
curves  such  as  shown  in  Figure  B-8. 

Dynamic  Analysis 

The  above  static  expressions  of  Equations  E-l  through  E-3  are  used  for  the  simplified 
dynamic  analysis  by  introducing  a  dynamic  amplification  factor,  D,  for  a  step  function  load, 
as  follows: 


D  =  1 


cos 


27ft 

T 


where  T  is  the  appropriate  response  period  of  each  of  the  three  functions  considered.  The 
amplification  factor,  D,  reflects  the  time-variant  characteristics  and  the  frequency  of  each 
function.  The  above  expression  can  be  expressed  in  terms  of  the  functions  of  a  variable, 
x,  as  follows: 


D(x) 


cos 


2  7T  t 
T  (x) 


(E-6) 


where  x  is  the  particular  variable  of  Equations  E-l  through  E-3. 

As  shown  in  Appendix  C,  the  frequency  of  the  axial  force  corresponds  to  the  "breathing" 
mode  of  a  complete  ring.  The  simplified  dynamic  response  of  the  axial  force  is  given  by  the 
expression 


Pt(t) 


[* 


cos 


2nt 

T(Pt) 


(E-7) 
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where  the  period,  T  (Pt),  is  given  by 


T(Pt)  =  2  n 


The  moment  due  to  the  symmetrical  loading  varies  at  the  frequency  of  the  first 
symmetrical  mode  of  vibration.  The  time  variant  response  of  moment  at  the  quarter  points 
can  be  expressed  as 


Mj(t) 


1 


- cos  TOjl]  Mi 


(E-8) 


where  the  corresponding  period  is  given  approximately  by 


T  (Mj) 


where  Am  is  given  in  Equation  E-5. 

The  moment  due  10  the  antisymmetrical  loading  varies  at  the  fiequency  of  the  first 
antisymmetrical  mode  of  vibration.  The  time-variant  response  cv£  moment  at  the  quarter 
point  under  the  inward  loading  can  be  expressed  as 


m2  (t) 


r,  2wt  ~ 

=  j^l  -  cos  T  (M2j 


M- 


(E-9) 


and  the  corresponding  period  is  given  approximately  by 


T(M2) 


where 


C 


-L_LJL 

-  i 


and 


n 


131 


The  total  time-variant  moment  is  obtained  by  summing  the  results  of  Equations  E-8 
and  E-9  at  each  time.  The  stresses  and  strains  at  the  quarter  point  under  the  inward  loading 
can  be  obtained  by  the  same  procedure  used  in  the  static  analysis. 

It  should  be  noted  that  the  maximum  dynamic  amplification  factor,  D(x),  of  Equation 
E-8  is  two  for  a  step  function  load. 

Presentation  of  Numerical  Results 

A  computer  program  {Program  E-l)  based  on  the  procedure  described  in  the  previous 
sections  was  used  to  obtain  the  numerical  solution  of  the  simplified  analyses.  Comparisons 
of  the  natural  periods  of  vibration  for  the  first  antisymmetrical  mode,  the  first  symmetrical 
mode,  and  the  second  symmetrical  mode  are  nude  between  the  simplified  analysis  and  the 
corresponding  values  of  Appendix  D,  and  are  presented  in  Table  E-l.  Close  agreement  is 
obtained  in  all  cases. 

Comparisons  between  the  simplified  analysis  and  the  more  rigorous  analysis  of  the 
dynamic  response  of  bending  moment  and  axial  force  as  given  in  Appendix  C  are  presented 
in  Table  E-2.  Similar  comparisons  for  static  values  are  presented  in  Table  E-3.  Compari¬ 
sons  of  the  dynamic  response  of  axial  forces  show  excellent  agreement.  The  difference  in 
the  static  values  of  axial  forces  is  a  result  of  the  elastic  support  conditions  introduced  in 
Appendix  B.  The  static  bending  moments  of  the  simplified  analysis  are  slightly  higher  than 
the  corresponding  results  of  Appendix  B.  T.'.ie  large  difference  (about  20  percent)  in  the 
dynamic  moments  is  due  to  the  damping  effects  of  the  system  considered  in  Appendix  C;  the 
simplified  analysis  did  not  consider  damping. 


Table  E-l.  Comparison  of  Natural  Periods  of  the  Simplified  Analysis 
to  the  Results  of  Appendix  D 


Load 
(lb/in. ) 

Antisymmetrical 

(msec) 

Symmetrical 

(msec) 

First  Mode 

First  Mode 

Second  Mode 

1 

2 

1 

2 

1 

2 

200 

17.0 

16.7 

6.3 

7.1 

mm 

3.3 

400 

17.3 

16.7 

6.4 

7. 1 

■H 

3.3 

600 

17.6 

16.7 

6.5 

7.1 

mm 

3.3 

1  -  Results  of  the  simplified  analysis 

2  -  Results  of  Appendix  D 
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Table  E-2.  Comparison  of  Maximum  Dynamic  Moments  and  Axial  Forces 
of  the  Simplified  Analysis  to  the  Results  of  Appendix  C 
(at  3/4 -point  of  the  arch) 


Load 
(lb/in. ) 

Simplified  Analysis 

Appendix  C 

Moment 
(in.  -lb) 

Time 

(msec) 

Axial 

Force 

(lb) 

Time 

(msec) 

Moment 
(in.  -lb) 

Time 

(msec) 

Axial 

Force 

(lb) 

Time 

(msec) 

200 

204, 000 

9.0 

19, 100 

1.6 

161,000 

7.2 

19, 000 

mm 

400 

415,000 

9.2 

38, 100 

1.6 

327,  000 

7.2 

38, 100 

IB 

600 

634,000 

9.4 

57, 200 

1.6 

442, 000 

8.0 

57, 100 

1.8 

Table  E-3.  Comparison  of  Static  Moments  and  Axial  Forces 
of  the  Simplified  Analysis  to  the  Results  of 
Appendix  B 


i 

Load 
(lb/in. ) 

Moments 
(in.  -lb) 

Axial  Forces 
(lb) 

1 

2 

1 

2 

200 

103, 500 

95,300 

9,600 

8,500 

400 

210,700 

194,200 

19, 200 

17,000 

500 

266, 300 

245,200 

24, 000 

21,300 

1  -  Results  of  the  simplified  analysis 

2  -  Results  of  Appendix  B 


Program  E-l.  Simplified  Analysis  of  Circular  Arches 


C  SIMPLIFIED  ANALYSIS  OF  CIRCULAR  ARCHES 

DIMENSION  T(3)«X(3)*FT(3)»E<4) 

1  FORMAT  ( I4.6E12.6.14) 

2  FORMAT  ( / 1 4  *  2 31  i  HALF  LOADED  ARCH  AT.F8.2.6H  LR/IN/I 

3  FORMAT ( /6H  T IME6X4HAX  F7X6HMI MAX ) 7X3HS T (T> 7X5HS I ( B » 8X4HS ( T I 8X 
1  4HSI31/) 

4  FORMAT  IF8«4,6E12.6I 

5  FORMAT  {/9X2HEA1UX2HMU1UX3HEAY9X3HSUM9X2HEI/! 

6  FORMAH/5H  K*  «c  12»6*jHPHIO«  *E12.6*3H  PCR*  *E  12 .6 . 5H  FAX«,E12.6i 

7  FORMAT  (SH  N*.E12.6.SH  KS=.E12.6.5H  MA«»E12.6l 

8  FORMAT  U>H  TN»*E12»6 *5H  TMS« »E12 . 6 . 5H  TMA«,E12.6) 

r  READ  CONTROL  CARD 

REAO  1.  NI,R,PHIO,PO.OP.TO*DT,NJ 
EI-O. 

EA»C. 

EAY«0. 

UH«0. 

PUNCH  5 
DO  11  L*1.4 

READ  1»  J.H,S,Y«E(LI .0 
A»H«3 

UM«UM4-A*490./(  1728  ••386.4) 

AY«A*Y*EIL) 

£I«EI+AY*Y*A*H*H*E<L»/12. 

EA«EA+A*E(L) 

EAY-FAY+AY 

11  PUNCH  li  L.EA.UM.AY.EAY.EI 
Cl-EAY/EA 
C2-D-C1 
E  I  =E  I -EA*C  1*0 
PUNCH  1,  L»£A,UM,C1»C2»EI 
P  I  =3*  141S927 
DO  32  I  =  1  * N I 
P=PO+OP 
X(1  I  «=P*R 
C=2.*PI/PHIO 
PCR=EI*(C*C-1.)/R**3 
FAX*1./(1.-P/PCR) 

S=FAX*(R*PHIO)**2 

X( 2 l =S*P/< (l.-(2./(3.*C)l**2)*216.) 

X  (  3  »  *S*P/  M  1  .-1./  IC*C  I)  *-32.  J 
T ( 1  l=2.*PI*R*SORT  CUM/EA  > 

5* S* SORT (UM/E I ) 

T I2»*2.*S/(9.*PI ) 

0*CC*0»1.51  UC*C-1.  1 
T I3»=CI*S/(2.»PI > 


6s 

RsPHJC'sPCRsFAX 

PUNCH 

7. 

X( 1 )*X(2 ) «X( 3) 

PUNCH 

8. 

T(ll.T(2)»TI3) 

Continued 
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Program  E-l.  Simplified  Analysis  of  Circular  Arches  (Contd) 


PO*P 
P«2,*P  ‘ 

PUNCH  2,  |  ,P 
PUNCH  3 
TOO. 

DO  31  J ■ 1 *N J 
XT  *TO«-OT 
00  21  K  *  1  *  3 

AF«1«-C0S(?»*P!*XT/T(K|  I 
21  FT(M«X<K)*AF 
XM*FT<  ? l+FT ( 3  I 
ST«FT( 1 1 /PA+XM«C1 /FI 
SB«FTU  I  /FA-XM«C2/FI 
STT  =  ST*E ( 2  ) 

STB=SB*F ( a | 

PUNCH  4.  Xt.M  ( 1  I  .XM.iT  .;5ctiTT  *STd 

31  TO*XT 

32  CONTINJF 
PAu.ir 
fnd 
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Appendix  F 

PHOTOGRAPHS  OF  TESTED  ARCHES 
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Figure  F-l.  Antisymmetrically  loaded  arch  AH-1  after  a  static  load  of  876  lb/ in. 
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Figure  F-3.  View  of  arch  AH-2  before  test. 


Figure  F-5.  View  of  arch  AH-3  before  test. 
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